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ABSTRACT 


A  user’s  guide  is  presented  for  a  computer  program  developed  to  aid  in  the 
design  of  sonic-fatigue-resistant  aircraft  structure.  The  program  employs 
matrix  methods  to  calculate  statistical  measurements  of  response  (deflection 
and  stress)  for  complex  structure  subjected  to  pressure  loads  random  in  both 
time  and  space.  The  program  is  in  two  phases.  Finite-element  methods  are 
used  in  the  first  phase  to  determine  structural  characteristics  such  as  flexibility, 
natural  frequencies,  and  modes  of  vibration.  In  the  second  phase,  a  cross¬ 
power  spectral  density  loading  function,  is  generated  and  combined  with  struc¬ 
tural  characteristics  to  compute  response.  Either  cross-power  spectral  density 
or  joint  statistical  moments,  including  second  spectral  moments  useful  in 
fatigue  analysis,  can  be  computed  for  response.  The  loading  function  models 
a  decayed  progressive  wave  typical  of  laboratory  noise  sources.  Different 
loading  functions  can  be  supplied  by  the  user,  because  the  program  is  constructed 
in  modular  form.  The  program  was  written  for  the  IBM  7094  computer  pri¬ 
marily  in  FORTRAN  IV  language  with  a  MAP  language  matrix  manipulation 
module. 
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I 

INTRODUCTION 

For  many  years  there  has  been  a  need  for  better  methods  of  sonic  fatigue 
analysis  for  aircraft  structure.  Current  methods  commonly  rely  on  available 
test  data  that  are  not  always  adaptable  to  new  designs.  A  design  tool  is  needed 
by  which  sonic  fatigue  performance  can  be  established  as  a  design  parameter. 

RANVIB  (Random  Vibration  Analysis  System  for  Complex  Structures)  is  a 
computer  program  developed  to  aid  in  sonic  fatigue  analyses.  The  program 
employs  matrix  methods  to  calculate  statistical  measurements  of  response 
(deflection  and  stress)  for  complex  structure  subjected  to  random  noise  fields. 

The  basis  of  the  analysis  method  is  a  combination  of  matrix  structural 
analysis  and  random-process  theory.  Load-displacement  relationships  are 
determined  by  finite-element  analyses,  and  vibrational  properties  are  deter¬ 
mined  by  computing  eigenvalues  and  eigenvectors.  These  characteristics 
are  used  to  form  structural  admittance  functions.  Force  loads  such  as  noise 
pressure  waves  are  based  on  using  cross-power  spectral  density  (cross  PSD)  on 
a  mathematical  model.  Response  is  related  to  the  loading  by  admittance  functions. 

The  program  is  written  primarily  in  FORTRAN  IV  language  for  the  IBM 
7094  computer.  However,  a  machine  language  module  is  used  for  matrix 
manipulations. 

The  computer  program  report  is  in  two  parts: 

1.  Part  I — Engineering  User’s  Guide 

2.  Part  II — Computer  Program  Description 

This  first  volume  is  a  guide  for  an  engineer's  use  of  RANVIB.  The  method 
of  analysis  used  in  the  program  is  described  in  section  II.  Section  III  describes 
the  use  of  the  program.  Sample  problems  for  application  to  the  program  are 
presented  in  section  IV. 

A  report  oriented  to  programmer  use  is  presented  separately  in  Part  11, 
(reference  1).  Development  of  the  theory  of  this  program  and  its  application  to 
specific  problems  are  presented  in  AFFDL-TR-68-44  (reference  2).  An 
earlier  study  that  used  portions  of  this  program  is  reported  in  reference  3. 
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n 

ANALYSIS  METHOD 


The  method  of  analysis  combines  matrix  structural  analysis  techniques  and 
random  process  theory  to  determine  structural  response  to  random  pressure 
loads.  The  computer  program  is  modular  with  each  module  performing  a 
specific  function  in  the  analysis. 

As  shown  in  figure  1,  the  program  is  organized  into  two  separate  phases. 
The  first  phase  contains  the  structural  generation  and  the  eigenvalue  and 
eigenvector  solution  modules.  Excitation  generation  and  response  solution 
are  in  the  second  phase.  Information  from  phase  I  is  transferred  to  phase  n 
by  means  of  tape. 

1.  SCOPE  OF  APPLICATIONS 

This  program  will  help  the  analyst  design  structure  that  is  resistant  to 
sonic  fatigue.  This  is  done  by  providing  information  on  the  random  stress 
response  distribution  on  small  sections  of  complex  structure  responding  to 
random  pressure  loads.  These  stress  distributions  can  be  used  to  determine 
where  regions  of  intense  stress  occur  and,  hence,  provide  stress  information 
necessary  for  fatigue  life  prediction.  Knowledge  of  the  stress  distribution 
allows  the  analyst  to  redesign  to  relieve  these  stresses  and,  thereby,  provide 
increased  fatigue  life.  Used  in  parallel  with  tests,  this  program  can  provide 
design  charts  on  complex  structural  configurations  and  minimize  the  need  for 
large  numbers  of  costly  tests. 

The  program  will  accommodate  large  complex  structural  systems  that 
exhibit  linear  elastic  behavior  and  are  represented  by  beam  and  plate  elements. 
Features  of  the  program  include  a  general  three-dimensional  beam  element, 
quadrilateral  and  triangular  plate  elements  (isotropic  or  orthotropic)  with  both 
in-plane  and  out-of-plane  stiffiiesses,  element  in -plane  preload  capability,  and 
elastic  support  capability.  Inertia  characteristics  are  determined  by  masses 
lumped  at  discrete  points . 
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Applied  loading  is  a  force  that  fluctuates  randomly  in  both  space  and  time. 
It  is  assumed  that  the  loading  is  a  stationary,  ergodic,  random  process.  The 
cross  PSD  of  forces  applied  at  structural  node  points  in  the  direction  of  nodal 
deflections  is  used  to  describe  the  excitation. 


PHASE  I 


PHASE  II 


Figure  1.  Computer  Program  Organization 
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The  excitation  subroutine  used  in  this  report  describes  a  mathematical 
model  of  a  laboratory  noise  source.  This  subroutine  can  be  replaced  with  one 
that  simulates  any  other  excitation  source.  The  present  form  of  this  sub¬ 
routine  restricts  use  of  this  program  to  rectangular  panels. 

Calculated  responses  include  deflections  and  stresses.  These  stresses 
are  actually  internal  forces  in  the  structural  elements,  and  element  geometry 
must  be  used  to  determine  true  stress  (force  per  unit  area).  Either  response 
covariance  or  cross  PSD  can  be  determined.  In  addition,  stress  second 
spectral  moments  can  be  used  to  find  the  expected  number  of  zero  crossings  in 
stress  time  history  for  use  in  predicting  fatigue  life.  Calculations  for  expected 
number  of  zero  crossings  are  valid  only  when  response -level  distributions  are 
Gaussian. 

Analysis  options  for  various  forms  of  system  damping  are  available  in  the 
program.  A  general  analysis  option  allows  damping  of  a  general  viscous  form. 
Simplified  analyses  can  be  performed  if  certain  assumptions  are  made  about 
damping. 

2.  PHASE  I— STRUCTURAL  AND  VIBRATIONAL  ANALYSIS 

Structural  stiffness,  flexibility,  and  stress-deflection  matrices  are 
generated  in  the  first  phase  of  the  program.  The  flexibility  matrix  and  mass 
data  are  then  used  to  determine  the  natural  frequencies  and  modes  of  vibration. 
The  results  of  the  first  phase  can  be  examined  to  determine  the  adequacy  of  the 
structural  representation  before  proceeding  to  the  response- solution  phase. 

a.  Structural  Matrix  Generation 
(1)  General  Information 

The  program  module  that  generates  stiffness,  flexibility,  and  stress 
matrices  was  adapted  from  a  larger  matrix  structural  analysis  program 
(reference  4).  This  program  is  capable  of  analyzing  large  complex  structural 
systems  by  the  stiffness  or  displacement  method. 

The  structure  is  described  in  Cartesian  coordinates  by  a  set  of  control 
points  (nodes)  connected  by  plates  and  beams.  Each  node  is  assigned  6  degrees 
of  freedom:  rotations  about  each  of  the  coordinate  axes  and  displacements  in 


the  directions  of  the  coordinate  axes.  The  constraint  conditions  for  each 
degree  of  freedom  are: 


(a)  Free 

(b)  Fixed  (zero  displacement) 

( c)  Sprung  (a  real  or  imaginary  spring  attached  to  a  ground  point) 

Either  quadrilateral  or  triangular  plates  may  be  used  for  structural 
representation,  but  best  stress  results  are  obtained  for  rectangular  elements. 
Beams  may  be  straight  or  curved  and  may  have  uniform  or  nonuniform  section 
properties . 

The  normal  assumptions  in  structural  analysis  are  made,  i.e. : 

(a)  The  material  is  perfectly  elastic 

( b)  The  deflections  are  sufficiently  small  compared  to  the  size  of  the  structure 
that  secondary  deflections,  caused  by  interaction  between  the  applied 
forces  and  primary  deflections,  are  negligible. 

This  program  is  limited  to  the  evaluation  of  structures  that  can  be 
adequately  described  by  no  more  than  2,  000  nodes.  The  total  number  of  plates 
and/or  beams  is  limited  only  by  this  restriction.  Further,  the  nodes  are 
grouped  into  partitions.  The  maximum  number  of  rows  of  partitions  is  200, 
and  the  maximum  number  of  nodes  per  partition  is  10.  The  total  number  of 
partitions  in  the  stiffness  matrix  is  limited  to  800. 

(2)  Fundamental  Matrix  Equations 

The  concepts  of  the  matrix  structural  analysis  are  contained  in  texts  such 
as  reference  5.  A  few  of  the  fundamentals  are  presented  herein. 

Most  structural  engineers  are  familiar  with  deflection  influence  coefficients 
a  that  relate  the  deflections  6  at  one  point  on  a  structure  to  external  forces 
F  applied  at  that  and  other  points.  For  any  point  on  a  structure,  an  equation 
may  be  written  in  the  form 
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If  all  n  equations  for  this  structure  are  written,  the  resulting  set  of 
equations  may  be  collected  and  expressed  in  matrix  form  as 
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This  matrix  of  deflection  influence  coefficients  is  called  the  flexibility  matrix. 


A  less  familiar  concept  is  that  of  the  force  influence  coefficients  K  that 
relate  the  external  force  at  one  point  on  a  structure  to  deflections  at  that  and 
other  points.  For  any  point  on  a  structure,  an  equation  may  be  written  in 
the  form 


F.  =  K  6  +  K  6  +  .  .  .  +  K.  .  6.  +  .  .  .  +  K.  6 

l  ll  1  i2  2  li  l  m  n 


(4) 


The  collection  of  these  force  equations  that  is  comparable  to  the  flexibility 
matrix  is 


or 


F1 

"Kll 

C\J 

*  *  *  Kln 

F2 

= 

K21 

K22  * 

•  • 

*  *  *  K2n 

•  •  •  • 

F 

n 

Knl 

4: 

.  .  .  K 

nn 

{F}  =  [K]  {6} 


(5) 


(6) 


The  matrix  of  force  influence  coefficients  is  conventionally  called  a 
stiffness  matrix,  and  the  individual  terms  in  the  matrix  are  called  the  stiffness 
coefficients.  The  stiffness  matrix  is  the  inverse  of  the  flexibility  matrix. 
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(3)  Matrix  Generation 


Generation  of  the  stiffness  matrix  is  accomplished  by  generating 
individual  element  matrices  and  merging  them  into  a  complete  structural 
stiffness  matrix.  Generation  of  element  stiffnesses  is  discussed  briefly  in 
appendix  I.  The  total  stiffness  matrix  is  modified  to  include  boundary  con¬ 
ditions  of  the  supported  nodes. 


Not  all  nodal  displacement  freedoms  are  retained  in  the  dynamic-response 
analysis  performed  by  the  RANVIB  program.  Thus,  the  structural  matrices 
must  be  compatible  with  the  retained  freedoms.  The  so-called  "reduced” 
stiffness  matrix  has  the  apparent  effect  of  selecting  a  smaller  set  of  nodal 
freedoms  in  the  matrix.  Reduced  freedoms,  which  are  assumed  to  be  unloaded 
by  either  applied  or  inertia  forces,  do  not  appear  explicitly  but  are  taken  into 
account.  The  complete  stiffness  matrix  is  rearranged  and  partitioned  so  that 
the  retained  freedoms  are  separated  from  the  unloaded  freedoms  to  be  reduced. 
Then  equation  (6)  can  be  written  as 


F1 

II 

i 

i 

^12 

o 

-  K21 

~  r 

i 

i 

K22 

The  reduced  displacements  are  given  by 


(?) 


Therefore, 


i&2)  =  -  [K 22]_1  [K21]  {61} 


{F 


!>  =  -  [K12]  [K22]  1  [K21]]  {6x} 


(8) 

(9) 


The  reduced  stiffness  matrix  is 


[Kr]  =  [K11]  “  [K12]  [K22]_1  [K21]  (10) 

The  corresponding  reduced  flexibility  matrix  is  equal  to  the  inverse  of  |  Kr  ] 

When  nodal  displacements  are  obtained,  they  can  be  used  to  obtain  internal 
forces  (stresses)  in  the  individual  members.  The  stress  deflection  relation¬ 
ships  can  be  written  in  matrix  form  as 


{a}  =  [S]  {6} 


(11) 
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where  { cr }  is  a  column  of  member  stresses  and  [S]  is  the  stress -deflection 
transformation  (stress  )  matrix.  The  stress  matrix  derived  in  the  program  is 
a  reduced  matrix  determined  in  the  same  manner  as  [Kr]  .  Equation  (11)  can  be 
written  as 


{a}  =  [S1  |  S2] 


(12) 


2 


Substitution  of  equation  (8)  into  equation  (12)  results  in 


{o}  =  [[Sl]  -  [S2]  [Kg2 1_1  tK21]]  {61}  (13) 

Therefore,  the  reduced  stress  matrix  is  given  by 

[Sr]  =  [S1]  -  [S2]  [K22]_1  [K21]  (14) 


For  the  remainder  of  this  report,  the  r  subscript  is  deleted  and  the 
structural  matrices  used  refer  to  reduced  matrices. 

Some  problems  require  an  analysis  that  will  include  the  effects  of  large 
axial  (or  in-plane)  loads.  Obviously,  these  can  cause  a  significant  change  in 
stiffness  of  the  structure.  This  effect  can  be  included  by  using  the  so-called 
geometric  stiffness  matrix.  Hence,  a  complete  stiffness  matrix  can  be 
formed  so  that 

[K]  =  [Kq]  +  [Kx]  (15) 


where:  [Kq]  =  the  usual  stiffness  matrix 

[  K11  =  the  geometric  stiffness  matrix 

Now  equation  (6)  can  be  written  by  substituting  equation  (15).  All  subsequent 
equations,  partitions,  etc. ,  are  valid  after  this  substitution.  For  a  more- 
complete  discussion  of  this  concept,  see  reference  6.  A  description  of  the 
[Ki]  matrices  is  contained  in  appendix  I. 

b.  Vibration  Analysis 

The  vibration  analysis  module  of  the  program  determines  natural 
frequencies  and  modes  of  vibration  for  the  system.  These  vibrational  char¬ 
acteristics  are  determined  from  the  undamped,  unforced  equation  of  motion: 

[K]  { 6(t ) )  +  [M]  {6(t)J  =  0  (16) 
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where  t  denotes  time.  The  mass  matrix  [M]  contains  inertia  coefficients 
corresponding  to  the  retained  freedoms  of  the  idealized  system. 


Solutions  of  equation  (16)  are  assumed  in  the  form 

(6(t)}  =  { 4> ^ ^  ^  1  exp  +  0)} 

J 

where:  {  0^}  =  a  vector  describing  a  deflected  shape 
co  j  =  a  circular  frequency 

ip  =  a  phase  angle 

Substitution  of  equation  (17)  into  equation  (16)  results  in  the  classical 
eigenvalue  problem 

|  [D]  -  A [I]  |  =  0 

where  [I]  is  the  identity  matrix. 

The  dynamic  matrix  [D]  is  given  by 

[D]  =  [K]-1  [M] 

and 

X  -  -i 
0) 


(17) 


(18) 


Each  eigenvalue  Aj  determines  a  natural  circular  frequency  coj  .  Associated 
with  each  eigenvalue  is  an  eigenvector  {00)}  that  represents  the  mode  shape 
for  that  frequency. 

The  eigenvalues  of  the  dynamic  matrix  are  solved  by  the  QR  method 
(reference  7)  in  this  program.  This  method  determines  all  eigenvalues  by  a 
series  of  transformations  of  the  dynamic  matrix.  In  general,  only  about  the 
lowest  25  percent  of  the  frequencies  and  modes  are  physically  realistic  because 
of  the  idealized  representation  of  the  system. 

The  flexibility  matrix  is  read  from  tape  in  this  module  and  the  mass 
matrix  is  read  from  cards.  To  provide  a  diagonal  mass  matrix,  the  inertia 
characteristics  are  represented  by  lumped  point  masses.  The  maximum 
matrix  size  is  of  100th  order,  and  up  to  25  modes  can  be  determined. 
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3.  PHASE  H — RANDOM  LOADS  AND  RESPONSES 


In  phase  n,  results  from  phase  I  and  information  obtained  from  the 
excitation  module  are  used  to  determine  responses  in  the  various  solution 
modules.  Phase  II  can  be  repeated  using  different  solution  options  without 
repeating  phase  I. 

a.  Random  Pressure  Loads 

(1)  General  Considerations 

The  RANVIB  analysis  system  was  constructed  for  use  on  structural 
systems  responding  to  a  class  of  fluctuating  pressures  called  sonic  loads.  The 
system  can  also  be  used  when  there  are  other  types  of  external  forces.  The 
only  requirements  are  that  the  applied  forces  be 

(a)  ergodic,  stationary,  random  processes,  and 

(b)  independent  of  structural  response  level  (aerodynamic  forces  excluded) . 

Matrices  describing  applied  forces  are  generated  by  a  random-loading 
module  (called  RANLOD)  in  the  phase  II  program.  This  loading  module 
describes  convective,  random,  fluctuating  pressures  by  use  of  a  mathematical 
model  based  on  decayed  progressive  waves.  This  form  of  loading  frequently 
occurs  in  sonic  test  laboratories  where  reverberation  on  reflection  is  small. 
The  output  of  RANLOD  is  force  cross-PSD  matrices.  RANLOD  can  be  replaced 
by  user-supplied  modules  that  describe  other  types  of  loadings  meeting  the 
general  requirements  stated  above.  The  loading  module  for  attached  turbulent 
boundary  layers  described  in  reference  3  is  an  example  of  the  use  of  another 
type  of  random-loading  module.  The  following  section  describes  the  loading 
module  for  the  decayed  progressive  wave;  it  serves  as  an  example  of  what  must 
be  considered  if  the  user  wishes  to  provide  his  own  RANLOD. 

(2)  Decayed  Progressive  Wave 

Sound  pressure  on  structural  panel  sections  located  at  great  distances 
from  any  type  of  random  sound  source  can  be  approximated  by  decaying  plane, 
progressive  waves.  An  example  of  a  progressive  wave  is  shown  in  figure  2. 
The  wave  front  will  generally  intersect  with  the  structural  surface  at  an  oblique 
angle  so  that  trace  wave  fronts  propagate  across  the  panel  forming  what  is 
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called  progressive  waves.  Normal  incidence  (plane  waves  parallel  to  the 
surface  of  a  flat  panel)  and  grazing  incidence  (plane  waves  perpendicular  to 
the  surface  of  a  flat  panel)  are  special  cases. 

Force  cross-PSD  matrices  [$F(ico)]  are  generated  and  used  in  the 
phase  II  solution  program.  These  matrices  are  Hermitian.  They  can  be 
described  as  the  sum  of  a  real  matrix  and  an  imaginary  matrix;  these  matrices 
are  force  co-power  [CF(co)]  and  quad-power  [QF(w)J  spectral  density 
(co-PSD  and  quad- PSD)  matrices: 

[4>F(iu)]  =  [Cp(aj)  ]  +  i[QF(w)]  (19) 

The  ico  denotes  that  cross  PSD  is  a  complex  function  of  frequency  w  . 


z 


Figure  2.  Pressure  Wave  and  Structural  Panel 
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The  sonic  type  of  distributed  pressure  loads  are  related  to  force  by 
the  following  equations: 


[Cjw)]  =  lA]  [C  U)]  fAj  (20) 

r  P 

[Qf(u3)  ]  =  CAJ  [Qp(u>)]  IX)  (21) 

where:  [  Cp(co) ]  =  pressure  co-PSD  matrix 

[Qp(cu)]  =  pressure  quad-PSD  matrix 

[  aJ  =  a  diagonal  matrix  of  areas  associated  with  each 
retained  node 

When  the  fluctuating  pressure  load  is  modeled  by  a  decayed  progressive 
wave,  the  terms  of  the  pressure  loading  matrices  are  (reference  2)  as  follows: 


C  (w)  =  <&(w)  exp  [-d  (€,  .2  +  ] 

P  .  .  1J  -L  J 

1J 

cos  [ — —  ( £ .  COS0  +  n..  sin0 ) ] 

ct  1J  1J 

2  2  1/2 

Q  (u>)  =  t’(w)  exp  [-d  +  )  ] 

Pij 

sin  [ —  (£.  .  cos0  +  n.  .  sin0 ) ] 

ct  1J  1J 

where:  $(cu)  =  pressure  PSD 
d  =  decay  constant 

t  =  x-direction  separation  distance  between  node 
points  i  and  j 

n  =  y-direction  separation  distance  between  node 

ij 

points  i  and  j 

0  =  propagation  angle  of  the  wave  front  over  the  panel 
Angle  9  and  trace  velocity  ct  are  shown  in  figure  3. 


(22) 


(23) 
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Figure  3.  Pressure-Wave  Propagation  on  Panel  Surface 


The  angle  0  defines  the  direction  that  the  trace  of  the  pressure  wave 
fronts  propagate  over  the  panel  surface.  This  angle  is  defined  by 


(24) 


where  cx  and  cy  are  the  respective  phase  velocities  of  the  pressure  wave 
in  the  x  and  y  directions  of  the  panel.  Wave-front  phase  velocities  in  the 
plane  of  the  panel  are  related  as  follows: 

c,  =  C  COS0 
t  X 

c.  =  c  sine 

t  y 

When  angle  0  approaches  zero,  a  more  accurate  definition  of  ct  results 

when  using  the  cosine  equation.  One  method  of  calculating  c  and  c  is 

from  knowledge  of  wave  front  transit  times  along  the  panel  edges  rx  and  r 

and  of  the  edge  lengths  a  and  b  .  Thus,  c  =  a/T  and  c  =  b/r 

x  x  v  '  v  ’ 


(25) 

(26) 
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The  loading  function  has  two  general  forms,  which  depend  upon  how  pressure 
wave  fronts  sweep  across  the  panel.  The  general  forms  are  as  follows: 

(a)  Normal-incidence  waves 

(b)  Progressive  waves 

Normal-incidence  waves  occur  when  wave  fronts  are  parallel  to  the  panel 
surface.  In  this  case,  the  trace  velocity  on  the  panel  approaches  infinity,  and 
the  trigonometric  functions  in  equations  (22)  and  (23)  have  constant  values: 

c°s  U..  cos8  +  n..  sine)]  =  i  (27) 


•  r  to)  / 

Sln  ( 5 •  .  cos  0  +  n .  .  sine ) ]  =  0 
-J-J  ■*  1 


1J 


(28) 


The  quad  PSD  is  everywhere  zero,  and  co-PSD  simplifies  to  the  product  of 
pressure  PSD  and  a  spatial  decay  term. 

A  progressive  wave  occurs  whenever  trace  velocity  c,  has  finite  value. 
There  are  three  types  of  progressive  waves: 

(a)  General  ( 0  t  c.  <  c  ,  c  <  c  ) 

d.  Xi  x  t  y 

(b)  Wave  progressive  in  the  x  direction  (0  =  nir,  c  =  c  ) 

t  x ' 

(c)  Wave  progressive  in  the  y  direction  ( 0  =  ^-2-n~IA7r  ,  c  =  c  ) 

2  ^  y 

where  n  is  any  integer.  When  the  pressure  wave  is  progressive  in  the  x  or 
y  direction,  respectively,  trigometric  terms  in  equations  (22)  and  (23), 
respectively,  reduce  to 


r  uj  ^  1 

cos  [—■£..] 
c 

X 

sin  [—  5.  ] 
Cx  1J 


(wave  progressive  in  the  x  direction) 


cos 

sin 


r  <*>  i 

r _ -1 


(wave  progressive  in  the  y  direction) 


15 


The  type  of  loading  function  generated  by  the  program  is  controlled  by  option 
and  parameter  inputs.  Selection  of  normal  incidence  or  progressive  wave 
forms  is  controlled  by  an  option  control  card.  In  additon,  the  type  of  progressive- 
wave  loading  function  that  is  generated  is  controlled  by  the  input  phase-velocity 
parameters.  When  the  wave  is  progressive  in  the  direction  of  one  of  the 
coordinate  axes,  phase  velocity  in  the  perpendicular  direction  is  infinite.  A 
zero  is  used  to  indicate  this  infinite  phase  velocity  in  one  direction.  The 
presence  or  absence  of  a  zero -coordinate  trace  velocity  then  controls  selection 
of  the  type  of  progressive  wave  loading. 

The  pressure  spectra  characteristics  are  determined  by  <|>(w)  .  These 
values  are  card  inputs.  A  pressure  PSD  must  be  the  input  for  each  frequency 
w  used  in  the  solution  computations  in  phase  n.  These  frequencies  are  either 
card  inputs  or  automatic  inputs  from  the  phase  I  analysis  of  natural  frequencies. 
The  source  of  frequency  data  depends  upon  whether  or  not  normal  modes  are 
used  in  the  solution  analysis  and  upon  whether  a  response  PSD  or  joint  moment 
is  calculated.  Refer  to  section  II  3.b.  for  a  discussion  of  response  solutions. 


Separation  distances  are  calculated  from  coordinates.  The  separation 
distances  between  node  points  i  and  j  are 


nij 


(29) 


The  locations  of  node  points  are  determined  by  the  idealization  used  in  phase  I 
to.  simulate  the  structure.  Most  structures  analyzed  for  response  to  sound  are 
panel  sections  that  can  be  idealized  by  locating  all  node  points  on  the  set  of 
intersections  formed  by  groups  of  lines  parallel  to  the  x  and  y  axes  (figure  4). 
RANLOD  requires  that  this  structural  idealization  be  used  here  so  that  separa¬ 
tion  and  areas  for  pairs  of  node  points  are  calculated  automatically;  hence  the 
present  form  of  the  loading  module  restricts  program  use  to  rectangular  panels. 
The  advantage  of  this  is  that  the  total  set  of  coordinates  need  not  be  used.  The 
only  inputs  that  are  necessary  are  the  sets  of  line-to-origin  distances.  The 
area  A  for  node  k  then  is  calculated  automatically  using 

\  =  f  Ui+1  -  xi-l>  lyj+l  -  yj-l>  (30) 
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Figure  4.  Structural  Nodes  on  Rectangular  Grid 


This  set  of  areas  is  grouped  to  form  diagonal  area  matrices  f  A  J  used  to 
premultiply  and  postmultiply  pressure  relationships  as  shown  in  equation  (21). 

Each  node  point  of  a  structure  has  six  potential  freedoms  of  motion,  but 
the  freedoms  of  most  node  points  are  not  retained  in  a  structural  analysis. 
Certain  node  points,  such  as  those  representing  a  boundary  of  a  clamped  plate, 
are  totally  constrained  and  have  no  freedoms.  Other  nodes  may  have  some  or 
no  constraints  but  are  not  directly  loaded  by  an  applied  or  inertial  force.  The 
freedoms  of  these  nodes  need  not  be  retained  in  the  structural  analysis, 
because  the  stiffness  matrix  of  retained  freedoms  includes  their  effects. 
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The  freedoms  that  are  retained  are  of  node  points  that  are  directly  loaded. 
The  applied  force  from  sonic  or  small-scale  fluctuating  pressures  acts  normal 
to  the  surface  of  an  aircraft  skin;  the  pressure-loaded  finite  elements  of  these 
skin  panels  have  only  the  translational  freedoms  normal  to  the  panel  surface 
(vertical).  Thus,  the  pressure  loads  that  must  be  included  in  a  force  cross- 
PSD  matrix  are  normal  forces  corresponding  to  the  node  points  with  retained 
vertical  translation  freedoms.  The  program  requires  that  inertial  loads 
correspond  to  these  same  retained  vertical  freedoms.  This  is  because  the 
automated  separation  generation  feature  precludes  generating  a  cross- PSD 
loading  function  with  zero  load  on  any  retained  freedom. 

The  Hermitian  properties  of  [  $F(icu)  ]  can  be  used  to  simplify  the  form 
of  the  co-  and  quad- PSD  matrices.  The  Hermitian  property  is 


[$F(iuO]  =  [4>F*(ia,)]T 

where  the  asterisk  denotes  the  complex  conjugate,  and  [  ]  denotes  the 
transpose  of  the  matrix.  Terms  of  the  matrices  have  the  following  properties: 


cf..k> 

11 

=  CF..(uk)  =  W 

JJ 

(31) 

CF..<V 

1J 

=  °F..(V 

(32) 

QF..(wk) 

n 

=  Qf  (o>k)  =  0 

jj 

(33) 

v/v 

1J 

3 

•H 

gP" 

1 

II 

(34) 

The  force  co-  and  quad-PSD  matrices  are  constructed  by  ordering  terms  to 
match  the  order  in  the  structural  stiffness  matrix  developed  in  phase  I.  The 
force  co-PSD  matrix  is 


0(wk) 

C12('V 

C13(“k> 

-  -  '  Cln 

C12(V 

$(u>k) 

C23(\J 

-  -  -  C2n 

n 

3 

Ik 

o 

C13<“k> 

1 

C23^k^ 

1 

4>(“k) 

C3n 

1 

_ciA> 

1 

C2nk> 

1 

CJ\> 

- $(w 

(35) 
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The  force  quad- PSD  matrix  is 


0 

w 

i 

3 

on 

a?-1 

-  -  vv 

[QFUk)J  = 

-w 

0 

-  -  w 

-Q13(“k) 

1 

-CW<V 

0 

1 

-  -  w 

1 

1 

"W 

1 

-w  - 

1 

-  0 

(36) 


where  subscript  a  is  the  number  of  retained  freedoms.  For  simplicity,  the 
F  subscript  is  dropped  in  the  terms  of  the  matrices. 

Diagonal  terms  of  the  co-PSD  matrix  are  the  collection  of  PSD’s  at  all 
node  points.  The  pressure  on  all  node  points  is  assumed  to  be  approximately 
constant,  so  that  all  diagonal  terms  are  equal.  If  N  frequencies  are  to  be 
used  in  the  response  solution,  there  will  be  N  co-PSD  and  N  quad- PSD 
matrices. 


When  phase  velocities  cx  and  c^  are  neither  known  nor  readily  obtain¬ 
able,  it  may  be  necessary  to  calculate  c^  and  9  directly  from  knowledge  of 
the  direction  of  propagation  of  the  free- space  wave  fronts.  This  can  be  done 
by  using  direction  cosines  to  define  a  normal  to  the  pressure  wave  front. 
Angle  y  is  generally  known  as  the  angle  of  incidence.  These  direction 
cosines  are  shown  in  figure  5.  The  panel  is  in  the  xy  plane.  The  relations 
are 


COS0  = 

cosa 

siny 

(37) 

c  =  c 

cosa  ) 

X 

cosg  I 

(38) 

c  =  c 

y 

c  =  ct. 

siny 

(39) 

where:  c 

e 

Phase  velocities 
input  of  cx  and 


=  acoustic  velocity  of  free-space  wave  front 
=  angles  used  to  form  the  direction  cosines  of  vector  <T 
=  angle  c^  makes  with  the  x  axis 

cx  and  Cy  can  be  calculated  by  using  equation  (26),  and  the 
c  can  then  be  used  in  the  program  to  calculate  c  and  9  . 

y  t 
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Figure  5.  Pressure-Wave  Group  and  Phase  Velocities 
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b.  Random  Response 


(1)  Cross-Power  Spectral  Density  (Cross  PSD) 

Response  time  history  cannot  be  determined  because  pressure  loads  on 
structure  vary  randomly  with  time.  However,  since  it  is  assumed  that  the 
loads  are  ergodic,  stationary,  random  processes  that  can  be  described 
statistically,  certain  statistical  information  about  the  random  responses  can 
be  obtained  from  the  loading  description  and  structural  characteristics. 

One  quantity  useful  for  describing  random  responses  is  cross  PSD. 

Reference  2  shows  that  deflection  cross  PSD  [  $  (iw)]  is  related  to  cross 

<5 

PSD  for  forces  acting  at  structural  node  points  [$F(iu>)]  by  the  expression 

[4>g(iu»)]  =  [H*(  iw)  ]  [$F(iu>)]  [H(iu>)]T  (40) 


where  the  diagonal  elements  of  [  4>^(iu;)]  are  the  spectral  densities  of  the 
deflections,  and  structural  characteristics  are  described  by  the  admittance 
matrix  of  complex  frequency-response  functions  (  H(iu>)]  .  The  admittance 
matrix  (reference  2)  is  given  by 

[H(iui)  ]  =  [-  w2  [M]  +  iu  [C]  +  [K]  ]_1  . 


Matrix  [C]  is  the  matrix  of  viscous  damping  coefficients.  Stress  response, 
or  internal  element  force,  is  related  to  deflection  response  by  stress-deflection 
relationships.  The  cross-PSD  matrix  of  element  stresses  [$  (icu)]  is 
given  as 


[<t>o(  ioo )  J  =  [sj  [*fi(iu)]  [S]T 


(42) 


Element  S..  of  matrix  [S]  is  the 
ij 

displacement  at  j  only. 


.th 

l 


elemental  stress  induced  by  a  unit 


(2)  Joint  Moments 

Other  information  on  random  responses  is  obtained  by  integrating  response 
cross  PSD.  The  integrals  of  cross  PSD  over  frequency  with  powers  of  fre¬ 
quency  as  weighting  functions  are  known  as  joint  spectral  moments.  When  the 
power  of  frequency  is  zero,  the  integrals  are  known  simply  as  joint  moments. 
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The  integral  of  a  deflection  PSD  gives  the  variance  or  mean- square  value  of 
that  deflection: 


(u)  do) 

Jo  6 


(43) 


The  rms  value  of  the  deflection  is  obtained  from  the  square  root  of  its 
variance.  Similarly,  an  average  value  can  be  obtained  for  the  products  of  the 
deflections  at  points  q  and  r  from  the  corresponding  cross-PSD  term: 


6  6 
q  r 


<J>  ( iio )  da) 


(44) 


All  such  terms  can  be  considered  simultaneously  and  written  as  a  matrix 
integration: 


[66] 


(io>)]  dw 


(45) 


The  elements  of  the  joint- moment  matrix  in  equation  (45)  are  space  cross¬ 
correlations  or  covariances  of  the  deflection  responses.  The  diagonal  elements 
are  the  mean-square  deflection  values. 

The  stress  covariance  matrix  can  be  obtained  by  substituting  stresses  for 
deflections  in  equation  (45).  The  result  is 


[oo]  =  [$  (id))]  dw 

'  o 

(46) 

=  [S]  [66]  [S]T 


Square  roots  of  the  diagonal  terms  of  the  stress  covariance  matrix  give 
rms  stresses. 

Information  useful  in  a  fatigue  analysis  can  be  obtained  from  second  spectral 
moments.  For  a  stationary  Gaussian  random  process,  the  expected  number  of 
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zero  crossings  Nq  of  a  stress  time  history  is  related  to  the  second-  and 
zero-order  moments  by  the  expression 


N 

o 


rOO  — 

(D2  $  ((d)  d(D 

0  ° 

f-oo 

$  (cd)  d(D 

o  °  J 

1/2 


(47) 


The  stress  second  moments  are  obtained  from  the  matrix  equation 


[$  (iw)]  dw 
a 


[S] 


t$.(ia))]  dw 
o 


tSj' 


(48) 


The  number  of  zero  crossings  for  a  particular  stress  component  is 
obtained  by  substituting  corresponding  diagonal  elements  of  the  matrices 
given  by  equations  (46)  and  (48)  into  equation  (47). 

(3)  Analysis  Options 

The  program  has  three  basic  options  for  calculating  the  responses 
described  in  paragraphs  n  2.b.  (1)  and  n  2.b.  (2). 

(a)  General  Viscous  Damping  (Option  1) 

Option  1  is  a  general  analysis  in  which  it  is  assumed  that  damping  is 
viscous  and  the  fluctuating  pressure  has  an  arbitrary  frequency  distribution. 

In  solving  the  cross  PSD,  the  matrix  of  damping  coefficients  is  read  from 
cards,  and  the  admittance  matrix  is  formed  for  each  frequency  specified. 

Since  a  matrix  interpretive  module  is  used  in  the  matrix  calculations,  com¬ 
plex  quantities  are  expressed  by  their  real  and  imaginary  parts.  Therefore, 
real  numbers  only  are  used  in  matrix  operations.  Equation  (41)  then  becomes 

[H(iu>)]  =  [J(u0]  -  i  [l(<d)  ]  (49) 

where 

[L(  id)  ]  =  [[(DC]  +  [K-(D2  M]  [idC]"1  [K— a)2  M]]-1 

[  J(id)  ]  =  [wC]-1  [K-id2  M]  [L((d)  ] 
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The  deflection  cross  PSD  for  each  of  the  frequencies  is  found  in  terms  of 
co-PSD  [C^]  and  quad  PSD  [Q^]  from  equation  (40)  using  the  symmetry 
properties  of  the  admittance  matrix 

[*6(iu»)  ]  =  [C6(u))]  +  i  L  Q6  ( m )  ]  (50) 

where 

[C6]  =  [J]  [CF]  [J]  +  [L]  [CF]  [L]  +  [J]  [Qp]  [L]  -  [L]  [Qf]  [j] 
[Q6]  =  [L]  [Cp]  [J]  -  [J]  [Cp]  [L]  +  [J]  [Qp]  [J]  +  [L]  [Qp]  [L] 

If  stress  cross  PSD  is  desired,  equation  (42)  is  used  to  obtain 

[*(!«*>)]  =  [C  (u)]  +  i[Q  (a.)]  (51) 

°  a  a 

where 

[Ca(u>)]  =  [S]  [C6(co)]  [S]T 

[Qa(“)]  =  [s]  [Qfi(a.)]  [s]t 

If  joint  moments  are  to  be  computed,  the  cross-PSD  matrices  are  numerically 
integrated  using  the  trapezoidal  rule.  The  deflection  covariance  matrix  is 
found  approximately  by  the  expression 

[66]  =  [$g(iu)x)] 


N-l  U)  .  .  -<D  .  , 

+  I  1+1  J-l 

j=2 


[*fi(i“j)]  +  ""Z"1'  [$6(ia)N)] 


(52) 


where  N  is  the  number  of  frequencies  for  which  cross  PSD  has  been 
determined. 


To  obtain  reasonable  accuracy  for  covariance,  cross  PSD  must  be  defined 
at  an  adequate  number  of  frequencies  over  the  significant  response  frequency 
range.  Second-deflection  spectral  moments  can  be  determined  in  a  similar 
manner  by  multiplying  each  term  in  the  summation  on  the  right-hand  side  of 
equation  (52)  by  the  square  of  the  corresponding  frequency.  Stress  covariance 
and  second- moment  matrices  can  then  be  determined  by  using  equations  (46) 
and  (48),  respectively. 
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(b)  Normal  Modes  (Option  2) 

Solving  for  deflection  cross  PSD  using  option  1  generally  involves  a 
considerable  number  of  matrix  inversions  and  multiplications  and,  thus,  may 
be  quite  time  consuming.  If  certain  assumptions  are  made  about  the  form  of 
damping  present  in  the  system,  the  motion  equations  can  be  uncoupled  when 
normal  mode  shapes  are  used  to  describe  displacements.  The  resulting 
simplified  analysis  is  option  2,  which  is  used  in  the  solution  module  of  the 
program. 

Damping  is  assumed  to  be  proportional  viscous  or  structural,  or  a  combi¬ 
nation  of  the  two  forms.  Viscous  damping  is  proportional  to  stiffness  and/or 
inertia.  The  damping  matrix  is  given  by 

[C]  =  p[M]  +  A[K]  (53) 

where  n  and  A  are  proportionality  factors.  For  structural  damping,  the 
matrix  [  C  ]  in  equation  (41)  is  replaced  by 

<4c]  =  g[K]  (54) 

where  g  is  the  structural  damping  coefficient.  It  is  shown  in  reference  2 
that  these  forms  of  damping  allow  the  equations  of  motion  to  be  uncoupled  so 
that  the  matrix  inversion  for  the  admittance  matrix  in  equation  (41)  reduces  to 
the  inversion  of  a  diagonal  matrix: 

[H(iw)]  =  [  <}>  ]  [Hj(iu))]  L<j>]  (55) 


where 


H .  ( ioj) 
J 


1 


2  2 

M.  (w.  —  to  +  i  2C  .  0)  . ) 

J  J  J  J 


The  term  u>.  is  the  j^1  natural  circular  frequency  of  the  system,  each  column 
3 

of  [  0]  is  a  normal  mode  shape  vector  {  0U/  }  ,  and  the  modal  damping 
factor  is  the  fraction  of  critical  damping  of  the  jth  mode.  The  term  M^ 
is  the  generalized  mass  defined  by 


M. 

J 


[M]  4(j)} 


(56) 
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The  modal  damping  factor  f  is  related  to  the  damping  parameters  by  the 
expression 


2c,.  =  —  +  X  a) . 

J  w  j 

J 


(57) 


A  similar  approximate  relationship  for  small  structural  damping  is 


«  ■  2cj 


(58) 


The  admittance  matrix  in  equation  (55)  can  be  expressed  in  terms  of  its  real 
and  imaginary  parts: 


[H  (iu>)]  =  [D]  -  i  [E] 

J 


(59) 


where  the  elements  of  the  diagonal  matrices  are 


2  2 
(D  .  -  0)^ 


M 


j  (o).2-(d2)2  +  (2^.0). a))2 
J  J  J 


(60) 


Ej  M. 


2 £  .a>  .a) 
J  J 


j  (u  2-u>2)2  +  (2C.U).0))2 
J  J  J 


(61) 


Substituting  equation  (55)  into  equation  (40)  and  simplifying,  the  expression 
for  deflection  cross  PSD  becomes 

[  $  (iw)]  =  [C  (a))]  +  i[Q-(oi)]  = 

0  6  0 


m  m 

l  l 

j=i  k=i 


(DA  +  EiV  [cF(Jk)(a,)] 


( jk) 


( jk)  / 


+  D  \  ( [Qp  J  (<*>)]  +  [Qp  J  («*>)]  ) 


+  i  [\  E.  ([CF(jk)(o))]  -  [CF(jk)(a,)]T) 


+  (D  Dk  +  E  E^ )  [Q^^M]] 


(62) 
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where 


[CF(jk)U)]  =  U(j)}T  [cFU)]  u(k)}  {4>(k) }T 

[QF(jk)(03)]  =  (4>(J } }  {<|>(j)}T  [Qf(u))]  {<t>(k)}  {<)>(k)}T 

The  double  summation  in  equation  (62)  is  over  m  normal  modes,  where 
m  is  equal  to  or  less  than  the  number  of  degrees  of  freedom  used  to  represent 
the  system.  For  most  problems,  reasonable  accuracy  is  obtained  when  m 
is  substantially  less  than  the  number  of  degrees  of  freedom. 

For  lightly  damped  systems,  the  direct  modal  products  (j  =  k)  and  the 
cross  products  for  which  numbers  j  and  k  are  nearly  equal  are  the  most 
significant  contributors  to  the  solution  in  equation  (62).  Thus,  some  of  the 
modal  cross-product  terms  can  be  omitted  in  the  calculations  with  no  signif¬ 
icant  change  in  results  but  with  considerable  savings  in  computing  time. 
Specified  to  the  program  in  this  option  can  be  parameter  K  ,  which  identifies 
how  many  cross-product  terms  are  to  be  computed  for  each  of  the  modes. 

Then  only  terms  for  which  I  j  -k  |  <  K  are  considered  in  equation  (62). 

To  compute  deflection  covariance  in  this  option,  it  is  assumed  that 
response  occurs  primarily  in  the  regions  of  natural  frequencies  (narrow- 
band  response)  and  that  excitation  varies  slowly  near  these  frequencies  (broad¬ 
band  excitation).  In  this  case,  the  integral  in  equation  (45)  can  be  evaluated  by 
treating  the  excitation  cross  PSD  as  constant  compared  to  the  admittance 
terms.  Evaluation  of  the  "constant"  excitation  matrices  midway  between  the 
appropriate  natural  frequencies  and  integration  of  equation  (62)  results  in 


[66]  = 


mm/  ,  .,s  /■“ 

i  i  h°F  <vj  (Dj\  +  Ej\)d“ 

j=l  k=l  l  J  J0  J  J 


t  (l8F(Jk)(“jk)1  +  [VJk)(v)lT> 


+  1  [(‘CF(Jk)(Vl  -  [CF<Jk>'VlT) 


D  du 


D.  E.  du) 
k  j 


+  j 


(Dj  \  +  Ej  V  a“1 } 


(63) 
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where 


U) 


2 


The  integrals  in  equation  (63)  have  been  evaluated  for  small  damping  in 
reference  8,  and  results  are  summarized  in  appendix  II.  As  in  the  cross-PSD 
calculation,  only  a  specified  number  of  the  cross-product  terms  of  equation  (63) 
are  computed. 

Second  spectral  moments  of  deflection  can  be  obtained  in  a  manner  similar 

to  the  calculation  of  the  covariance  matrix.  The  only  difference  being  that 

2 

the  integrands  of  the  integrals  in  equation  (63)  are  multiplied  by  oj  .  These 
new  integrals  were  also  evaluated  in  reference  8  and  are  summarized  in 
appendix  II. 

Stress  cross-PSD  and  joint  moment  matrices  are  computed  in  this  option 
in  the  same  manner  as  in  the  general  analysis  (option  1).  The  corresponding 
deflection- response  matrices  are  premultiplied  by  the  stress  matrix  [S  ]  and 


T 


postmultiplied  by  [  S] 


(c)  Normal  Modes  With  Cross  Terms  Omitted  (Option  3) 


Option  3  of  the  response- solution  module  is  a  further  simplification  of 
the  normal-mode  analysis.  In  this  option,  all  modal  cross-product  terms  are 
omitted  in  the  calculation  of  response  cross  PSD  and  joint  moments.  The 
accuracy  of  this  analysis  is  usually  sufficient  for  very  lightly  damped  systems, 
and  computation  time  is  significantly  less  than  that  required  for  the  other 
analysis  options. 

Disregarding  cross-product  terms  in  equation  (62)  and  utilizing  the 
Hermitian  properties  of  the  excitation  cross-PSD  matrices,  the  equation  for 
deflection  cross  PSD  reduces  to  a  single  summation: 


l  (D  2  +  E.2)  [c  (0))]  {*(  J  )}T  (64) 

J  J  F 


Response  quad  PSD  cannot  be  obtained  in  this  simplified  analysis,  because 
the  resulting  imaginary  component  is  null. 
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Deflection  covariance  matrix  evaluation  is  also  simplified  when  cross- 
product  terms  are  ignored.  Omitting  cross-product  terms  from  equation  (63) 
and  evaluating  the  remaining  admittance  integral  results  in  deflection  covari¬ 
ance  being  expressed  as 

m 

I 


[Vr] 


j=l 


u(j)}  {^(j )  }T  [c  (u>  )]  u(j)}  — — — - 

J  4m  2W  3 


J 


J 


(65) 


The  resulting  imaginary  component  of  the  covariance  matrix  is  also  null  in 
this  analysis. 

An  expression  similar  to  equation  (65)  results  when  evaluating  deflection 
second  spectral  moments.  The  term  in  the  denominator  of  the  right-hand 
side  of  equation  (65)  is  replaced  by  u>j  for  each  term  of  the  summation. 

Stress-response  matrices  are  determined  from  the  corresponding 
deflection  matrices  in  the  same  way  as  that  in  the  other  analysis  options. 
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Ill 

USE  OF  PROGRAM 


1.  COMPUTER  CONFIGURATION  REQUIREMENTS 

The  RANVEB  programs  were  written  in  FORTRAN  IV  and  MAP  to  be  run 
under  the  7094  Mod  II  IBSYS  operating  system.  By  proper  input /output  defini¬ 
tions,  it  has  been  executed  successfully  under  a  coupled  7040/7094  system. 

2.  TIMING  AND  OUTPUT  ESTIMATES 

a.  Phase  I 

Execution  time  for  phase  I  of  the  program  depends  on  many  parameters , 
including  the  number  of  node  points,  total  number  of  kinematic  degrees  of  free¬ 
dom,  number  of  retained  freedoms,  and  number  of  structural  elements.  The 
exact  relationship  between  these  factors  and  computing  time  is  not  known.  How¬ 
ever,  based  on  previous  problems,  a  rough  time  estimate  can  be  obtained  from 
the  number  of  nodes.  For  problems  in  which  the  number  of  nodes  is  less  than 
about  150,  the  execution  time  in  minutes  is  equal  to  about  one  fourth  the  number 
of  nodes.  When  stress-deflection  matrices  are  not  computed,  the  time  is  about 
20-percent  less. 

For  larger  problems,  computing  time  is  greater  than  estimated  by  the  above 
linear  relationship.  However,  not  enough  information  is  available  to  establish 
a  guide  for  estimating  time  for  these  problems. 

The  amount  of  printed  output  also  varies  with  problem  size.  A  guide  for 
estimating  the  number  of  output  lines,  excluding  program  loading  maps,  has 
been  established.  Structural  input  data  are  printed  as  they  are  read,  creating 
about  one  line  of  output  for  each  node,  plate  element,  and  beam  element  in  the 
idealization.  Further  output  depends  on  the  number  of  retained  freedoms  n 
and  the  number  of  computed  mode  shapes  m  .  An  approximate  expression  for 
the  number  of  additional  output  lines  is 

50  +  n  [1  +  3  (f)  +  (g)]  +  m  [2  +  (|)  ]  (66) 

The  numbers  in  parentheses  should  be  rounded  up  to  the  next  whole  number. 
Optional  diagnostic  printing  of  elemental  structural  matrices  can  create  much 
additional  output. 
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b.  Phase  II 


Execution  time  for  phase  II  depends  upon  analysis  and  solution  options  as 
well  as  problem  size.  Most  of  the  information  available  is  for  the  option  3 
analysis.  For  deflection  and  stress  joint-moment  solutions,  an  approximate 
value  for  execution  time  is  0.  008  (m  n  )  minute.  The  time  is  about  30-percent 
less  when  second  spectral  moments  are  not  computed. 

Execution  time  (minutes)  required  for  an  option  3  cross- PSD  analysis  for 
N  frequencies  is  approximately 

0.004  nmN  (67) 

In  an  option  2  analysis,  cross-modal  response  terms  are  added  to  option  3 
responses  to  obtain  the  total  results.  Thus,  the  increase  in  execution  time 
over  a  corresponding  option  3  computing  time  depends  upon  the  number  of  cross 
products  included  in  the  analysis.  For  K  cross  products  for  each  mode, 
option  2  execution  time  is  roughly  (1  +  2K)  times  the  value  for  an  option  3 
analysis. 

Little  information  is  available  on  option  1  execution  times,  because  only 
small  problems  have  been  used  with  this  option.  The  time  for  a  cross-PSD 
solution  for  a  small  problem  was  about  the  same  as  an  option  2  analysis.  Com¬ 
putation  time  is  expected  to  be  quite  substantial  for  large  problems,  however, 
since  inversions  and  multiplications  of  n  by  n  matrices  are  performed  in  the 
analysis.  Also,  if  joint  moments  are  to  be  calculated,  deflection  cross  PSD 
must  be  computed  at  a  sufficient  number  of  frequencies  to  obtain  reasonable 
accuracy  in  the  numerical  integration  process. 

The  amount  of  printed  output  for  phase  II  also  depends  upon  the  analysis 
and  solution  options.  For  an  option  3  analysis  for  deflection  cross  PSD  at  N 
frequencies,  the  number  of  lines  of  output  data  is  approximately 

100  +  N[3  +  m  +  n  (y)  ]  (68) 

If  any  generated  force  cross-PSD  matrices  are  printed,  about  n(y)  additional 
lines  are  printed  for  each  matrix.  Stress  response  adds  15  lines  of  output  for 
each  beam  element  and  9  lines  for  each  plate  element  for  each  of  the  N 
frequencies  specified. 


32 


An  option  3  analysis  for  deflection  covariance  and  second  spectral  moments 
results  in  about  80  +  2n(^)  lines  of  output.  Stress  response  matrices  add 
about  30  lines  for  each  beam  element  and  18  lines  for  each  plate  element. 

Optional  printing  of  force  cross- PSD  matrices  again  adds  n(^)  lines  for  each 
matrix. 

Both  real  and  imaginary  parts  of  the  loading  and  response  matrices  are 
printed  for  an  option  2  analysis.  As  a  result,  an  option  2  analysis  creates 
about  twice  as  much  output  as  a  corresponding  option  3  analysis. 

The  number  of  lines  of  output  data  for  an  option  1  deflection  cross -PSD 
solution  at  N  frequencies  is  approximately 

80  +  3N  +  n  [2N  +  l]  <y)  (69) 

Optional  printing  of  force  cross  PSD  adds  2n(^)  lines  for  each  matrix.  If 
stress  cross  PSD  is  determined,  the  output  is  increased  by  30  lines  and  18  lines 
for  each  beam  and  plate  element,  respectively.  For  deflection  covariance  and 
the  second-spectral-moment  solution,  the  amount  of  output  is  the  same  as  that 
for  a  deflection  cross-PSD  solution  plus  4n(^)  lines.  Stress  covariance  and 
second-spectral-moment  responses  each  add  the  same  amount  of  output  as  that  of 
stress  cross  PSD. 

3.  INPUT  DATA 

a.  Control 

A  diagram  of  the  deck  setup  for  phase  I  and  phase  II  is  shown  in  figure  6. 
There  are  two  sets  of  control  cards  followed  by  the  problem  data  cards. 

The  first  set  of  control  data  is  related  to  the  computer  installation.  Cards 
such  as  job  number,  timing,  and  print  estimates  are  included  in  this  set.  Also, 
identification  of  tapes  for  the  program  and  input/output  must  be  specified. 

Control  cards  for  the  program  subroutines  that  are  on  tape  are  contained 
in  the  second  set  of  control  data.  Listings  of  these  control  cards  are  presented 
in  appendix  III,  pages  80  and  83  ,  for  phase  I  and  phase  II,  respectively. 
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Figure  6.  Deck  Setup 


b.  Input 

(1)  Phase  I  Data 

Input  to  phase  I  of  the  program  is  in  the  form  of  data  cards.  These  data 
describe  structural  information  necessary  to  generate  structural  matrices  and 
to  determine  vibrational  characteristics.  The  input  is  organized  into  the 
following: 

(a)  Control — structure  size,  partition  sizes,  and  output  options 

(b)  Nodal  data — nodal  coordinates  and  constraint  conditions  (fixed,  free,  etc.) 

(c)  Plate  data — node  numbers  and  properties  (moments  of  inertia,  areas,  etc.) 

(d)  Beam  data — node  numbers  and  section  properties 

(e)  Eigenvalue  data — matrix  size,  number  of  modes,  and  mass  matrix 

NOTE:  All  input  units  of  length  and  force  are  assumed  to  be  inches  and  pounds, 
respectively. 

There  are  two  types  of  coordinate  systems.  The  location  of  nodes  defining 
a  structural  assembly  are  defined  by  a  structural  coordinate  system  designated 
by  Cartesian  coordinates  with  lower-case  letters  x  ,  y  ,  and  z  .  Each  element 
has  its  own  local  coordinate  system  defined  by  cartesian  coordinates  designated 
by  capital  letters  X  ,  Y  ,  and  Z  .  These  local  coordinate  systems  are  shown 
in  figure  7  and  are  defined  as  follows: 
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NODES  3  AND  4  CAN 


(a)  Quadrilateral  Plate  fa)  Triangular  Plate 


— Jb  \+- 


SHEAR  AREAS:  Axz  =  tweb  .  d 
AXy  =  tf)  •  b  •  2  • 

(c)  Straight  Beam 


SHEAR  AREAS:  Axz  =  tf)  •  b  •  2  • 
AXY  =  lweb  ‘  d 


(d)  Curved  Beam 


Figure  7.  Local  Coordinate  Systems 


(a)  Plates — The  local  X  axis  originates  at  node  1  and  is  oriented  in  the 
direction  of  node  2.  The  local  Y  axis  originates  at  node  1  and  is  oriented 
perpendicular  to  the  local  X  axis  in  the  plane  of  nodes  1,  2,  and  3  in  such  a 
way  that  node  3  has  a  positive  Y  value.  The  local  Z  axis  is  established  by 
the  right-hand  rule. 

(b)  Straight  beams — The  local  X  axis  originates  at  node  1  and  is  oriented  in 
the  direction  of  node  2.  The  local  Z  axis  originates  at  node  1  and  is 
oriented  perpendicular  to  the  local  X  axis  in  the  plane  of  nodes  1,  2,  and 
3  in  such  a  way  that  node  3  has  a  positive  Z  value.  The  local  Y  axis  is 
established  by  the  right-hand  rule.  If  the  3rd  node  is  zero,  the  local  Z 
axis  is  assumed  to  originate  at  node  1  and  be  parallel  to  the  structural  z 
axis. 

(c)  Curved  beams — The  local  X  axis  originates  at  node  1  and  is  oriented  in 
the  direction  of  node  2  tangential  to  the  curve.  The  local  Y  axis  originates 
at  node  1  and  is  oriented  perpendicular  to  the  local  X  axis  in  a  radially 
outward  direction  in  the  plane  of  curvature  of  the  beam.  Node  3  defines 

the  plane  of  curvature  and  must,  therefore,  be  in  that  plane  and  can  be 
located  in  the  half  plane  (created  by  a  straight  line  connecting  nodes  1  and 
2)  not  containing  the  beam.  For  this  case,  enter  the  radius  of  the  curve 
as  positive.  Node  3  can  be  located  in  the  half  plane  containing  the  beam 
but  not  in  the  segment  between  the  beam  and  the  line  connecting  nodes  1 
and  2.  For  this  case,  enter  the  radius  of  the  curve  as  negative.  The  local 
Z  axis  is  established  by  the  right-hand  rule.  The  local  coordinate  system 
at  node  2  is  established  by  moving  the  coordinate  system  from  node  1  to 
node  2  along  the  curved  centerline  of  the  beam  and  maintaining  the  local 
Y  axis  in  a  radially  outward  direction. 

The  following  shows  the  input-data  preparation  procedures  for  phase  I: 
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Card  No.  Card  title  Column  Description  Format 

1  Option  control  card  01-10  Program  Option  I  10 
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1 — Print 

NOTE:  The  elemental  stress -matrix  print 

option  creates  much  additional  output 
and  is  normally  not  used. 
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2  (Cont. )  Spring  constants  13-24  Spring  constant  (in. -lb/rad)  for  rotation  6  E12.  4 
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(a)  Sandwich  Plate 


(b)  Stiffened  Plate 


Figure  8.  Sandwich  and  Stiffened  Plates 
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(Leave  blank  if  section  property  option  =  0  or  1. ) 
25-28  Print  option 
0 — No  print 

1 — Print  coordinate -transformation,  stiffness, 
and  stress  matrices  for  this  beam  element. 
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Geometric  stiffness-  01-12  Axial  load  in  the  local  X  direction  applied  on  E12.4 

matrix  input  (Omit  this  card  the  beam  (lb) (tension  positive) — Use  this  field 

if  card  2,  columns  47  and  48,  if  rotation  will  occur  about  the  local  Y  axis, 

are  either  blank  or  zero. ) 


2 A  (Cont. )  Geometric  stiffness-  13-24  Axial  loads  in  the  local  X  direction  applied  E12.4 

matrix  input  on  the  beam  (lb)  (tension  positive) — Use  this 

field  if  rotation  will  occur  about  the  local  Z  axis . 
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(See  figure  7.) 

25-36  A  area  carrying  axial  load  (sq  in.)  E12.  4 


4A  (Cont. )  Constant  cross  section  37-48  I  moment  of  inertia  about  Z  axis  (in.  )  E12.  4 


cm"  cm" 

tH  i— ( 

W  W 


'f  ^ 

I— I  t— I  I— I  I— I  ^ 


© 

u 

3> 


© 

0) 

m 


g’  £ 


D* 

CD 

(1) 

Cl 

cd 


N 

X 


Ci 

0) 

Id 

a 

cd 

Ci 

cd 

a 

CD 

c 

cd 


cd  .2 

m 

<D 

Cl 

Cl  rl 

«d  3 

CD  g 

•G  .2 

W  m 

Ci 

N  O 


O  CM 
CO  t- 

I  I 

CT>  r-l 

^  CD 


CD 

Cl 

cd 


o 

Cl 

cl 

cd 

o 

© 

> 


© 

o 

Oi 

Ci 

o 

o 

Ci 

© 

bO 

© 

X 

G1 

a 

m 

© 

•H 

tJ 

© 

a 

o 

Ci 

a 

G 

O 

•H 

4-> 

O 

0 

m 


u 

d 

0 

0 

0 

o 

c 

b£) 

•»-« 

0 

£ 

ctf 

0 

a 

d 

0 

c3 

a 

© 

m 

o 

Ci 

© 

N 


© 

> 

cd 

be 

© 

G 

© 

XI 


G 

O 

3 

a 

Ci 

o 

•a 


Ci 

© 

X 

a 

3 

G 


Ci 

© 

a 

o 

Ci 

Pi 


■d* 

o 


J*  < 

I  I  I 

rH  CM  CO 


N 

N  J* 

hh  <; 


LO  CD 


£ 

£ 

£ 

£ 

o 

o 

1—4 

o 

f— H 

o 

*0 

0 

0 

"© 

42 

42 

42 

X 

c3 

ci 

d 

cd 

i—4 

3 

3 

3 

3 

3 

3 

u 

a 

c, 

3 

a 

d 

i—4 

o 

o 

o 

o 

*4-4 

*4-1 

*4—1 

*4-1 

42 

0 

0 

0 

0 

0 

0 

0 

0 

0 

w. 

C/D 

CO, 

d 

0 

W 

vA 

s 

£ 

00 

CM 

CO 

o 

o 

1 

t-H 

t-H 

CM 

• 

CM 

1 

in 

Oi 

00 

i 

t> 

r-H 

o 

o 

i— H 

rH 

CM 

'O 

0 

3 

r— 1 

3 
0 
I  i 

d 

o 

0 

0 

•H 

t: 

© 

a 

o 

Ci 

a 

G 

O 

•H 

o 

0 

GO 


{3 

O 

•H 

o 

© 

.CO. 

_cd 

"3 

a 

Ci 

o 


bo 

G 

•H 

X5 

U 

O 

O 

o 

oi 


C 

o 

•H 

& 


Ci 

© 

a 

o 

Ci 

a 


3 


48 


25-36  A  (See  formula  below. )  E12.  4 

37-48  B  (See  formula  below. )  E12.  4 
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(2)  Phase  II  Data 


Phase  II  input  is  from  both  tape  and  cards.  The  phase  I  output  tape  must  be 
mounted  on  logical  unit  10  for  use  in  phase  II.  Additional  information,  which 
comes  from  data  cards,  is  organized  into  two  sections:  control  and  excitation 
descriptions. 

The  first  input  section  (control)  contains  analysis  and  response  solution 
options  and  problem-size  parameters.  Damping  values  and  frequency  informa¬ 
tion  are  also  included  in  this  section. 

The  second  input  section  (excitation  description)  contains  data  necessary 
for  the  generation  of  force  cross-PSD  matrices  for  a  panel.  Care  must  be  taken 
so  that  the  order  of  terms  in  these  matrices  duplicates  the  order  in  the  structural 
matrices.  Compatibility  can  be  obtained  by  taking  the  following  steps: 

(a)  Position  the  panel  in  the  first  quadrant  of  a  rectangular  coordinate  system 
such  that  the  first  structural  node  with  retained  deflection  is  nearest  the 
origin  (figure  9).  Indicate  by  vector  c^  the  direction  of  progressive  wave 
propagation  across  the  panel  relative  to  this  orientation. 

(b)  Determine  c  and  c  ,  the  x-  and  y-direction  phase  velocities.  Signs  of 

x  y 


(c)  Determine  the  direction  (x  or  y)  in  which  the  structural  nodes  are  ordered 
by  an  increasing  node  number.  This  direction  specifies  the  cyclic  direction 
of  node  numbering  (figure  11). 

(d)  Assign  an  origin-to-line  distance  for  each  line  of  nodes  (figure  12).  Dis¬ 
tances  for  two  lines  defining  panel  boundaries  for  each  direction  must  be 
included. 

Pressure  PSD  must  be  read  in  for  each  frequency  at  which  force  cross  PSD 
is  to  be  generated.  These  frequencies  depend  on  the  analysis  and  response 
options  being  used.  For  analysis  option  1  and  response  cross-PSD  solutions, 
pressure  PSD’s  correspond  to  the  frequencies  specified  in  the  first  input 
section.  When  joint  response  moments  are  calculated  in  analysis  option  3, 
pressure  PSD’s  must  be  specified  for  the  system’s  natural  frequencies  used  in 
the  analysis .  For  an  option  2  joint -moment  calculation,  PSD’s  are  read  in  two 
sets .  The  first  set  is  for  the  natural  frequencies  and  the  second  set  is  for  the 
cross -modal  terms.  These  cross -modal -term  frequencies  are  in  the  order 
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i— 1,2,  .  .  .  ,  m— 1 

j  =  i+1,  .  .  .  ,  i+K<m 

where  K  is  an  input  parameter  defining  the  number  of  cross  terms  used  in  the 
calculation  and  m  is  the  number  of  natural  frequencies  used  in  the  analysis. 
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CYCLIC  ORDERING  OF  RETAINED  NODES 
TYPICAL  NODE  LINES 


NOTE:  IN  THIS  EXAMPLE, 
THE  ORDERING  OF 
NODES  IS  CYCLIC  IN 
THE  y  DIRECTION. 


Figure  11.  Cyclic  Ordering  of  Nodes 
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DISTANCES  IN  NONCYCLIC  DIRECTION 
-n2  NODE  LINES 


Figure  12.  Node  Line  Distances 


The  following  shows  the  input-data  preparation  procedure  for  phase  II: 
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(a)  Control 

Card  No.  Card  title  Column  Description  Format 
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Card  No.  Card  title  Column  Description  Format 

2  (Cont. )  Parameter  card  11-20  Mass  proportional  damping  factor;  see  F10. 0 

equation  (53). 

21-30  Structural  damping  coefficient;  see  equation  (54).  F10.  0 

31-40  Number  of  retained  freedoms — 90  maximum  I  10 
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(b)  Excitation  Description 
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in  y  direction  (Leave  blank  when  this  phase 
velocity  is  infinite. ) 


Card  No.  Card  title  Column  Description  Format 

4  Grid-size  card  01-10  Number  of  coordinates  (including  two  boundaries)  I  10 

in  the  direction  of  cyclic  node  numbering 
11-20  Number  of  coordinates  (including  two  boundaries)  I  10 
in  the  direction  perpendicular  to  the  cyclic  node 
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4.  OUTPUT  DATA 


a.  Standard 

Input  data  for  the  structural  matrix  generation  module  in  phase  I  are 
printed  as  they  are  read.  The  generated  reduced  stiffness  matrix  is  printed 
only  if  the  print  option  is  exercised. 

Flexibility,  dynamic,  and  input  mass  matrices  are  printed  in  the  eigen¬ 
value  module.  After  eigenvalues  and  eigevectors  have  been  determined, 
matrix  {  [0]  T  [m][0]}  is  printed.  Then,  diagonal  elements  of  this  matrix 
(i.  e.  generalized  masses)  are  printed  separately.  The  final  output  of  the  first 
phase  is  the  mode- shape  vectors  and  natural  frequencies. 

Ordering  of  terms  in  the  output  matrices  corresponds  to  structural  node 
numbering,  whereby  the  first  element  corresponds  to  the  retained  freedom  with  the 
lowest  node  number.  If  more  than  one  freedom  is  retained  for  a  node,  displace¬ 
ments  for  that  node  are  in  the  order  (0V  ,  Q  ,0,6,  6  ,  6  )• 

v  X  ’  y  1  z  x  y  z' 

In  phase  n,  input  data  specifying  solution  options  and  parameters  are 
printed  as  they  are  read.  Also,  natural  frequencies  and  generalized  masses 
are  printed  as  they  are  read  from  the  phase  I  tape.  Printing  of  generated 
force  cross-PSD  matrices  is  optional. 

Further  output  is  dependent  upon  the  solution  and  analysis  options  selected. 
Either  deflection  cross-PSD  matrices  are  printed  for  each  frequency  specified 
or  the  deflection  covariance  matrix  is  the  output.  For  analysis  options  1  and  2, 
real  and  imaginary  parts  are  printed  for  these  matrices,  whereas  only  real 
parts  are  printed  in  option  3.  Deflection  response  matrices  are  of  order 
n  by  n,  whereby  n  is  the  number  of  freedoms  retained  in  the  analysis. 


Ordering  of  deflections  in  the  matrices  again  corresponds  to  the  structural 
node  numbering. 

Stress  cross-PSD  or  joint  moment  matrices  are  printed  if  the  stress 
response  option  is  exercised.  Stress  response  matrices  are  printed  for  individ¬ 
ual  structural  members.  For  each  plate,  an  8-by-8  matrix  is  printed  with 

elements  corresponding  to  stresses  (Fx,  Fy,  NxY’  ^X’  ^Y’  MX’  MY’  MXY^‘ 
These  stress  components  (figure  13)  refer  to  forces  per  linear  inch  at  a  section 
through  the  center  of  gravity  of  the  plate  element.  Two  6-by-6  matrices  are 
printed  for  each  beam  element,  one  each  for  the  stresses  at  each  end  of  the 
beam.  Beam  stresses  (figure  14)  correspond  to  end  forces  (Mx>  My,  Mz>  Fx, 

FY’  fz>- 

b.  Diagnostics 

Error  messages  possible  in  the  structural  generation  module  are  as 
follows: 

( 1)  NODE  XX  OUT  OF  SEQUENCE  — Indicates  bad  nodal  data;  check  data  deck. 

(2)  INCORRECT  NODE  NUMBER,  PLATE  NO.  XX — Indicates  plate's  node 
number  is  greater  than  the  number  of  nodes  in  the  structure;  check  data  deck. 

(3)  NEGATIVE  MODULUS  OF  ELASTICITY,  PLATE  NO.  XX— Indicates 
negative  Young’s  modulus  for  plate;  program  converts  value  to  positive. 

(4)  NEGATIVE  POISSON  RATIO,  PLATE  NO.  XX — Indicates  negative  Poisson 
ratio  for  plate;  program  converts  value  to  positive. 

(5)  INPUT  NOT  IN  PROPER  SEQUENCE,  PLATE  NO.  XX  SHOULD  BE  XX— 
Indicates  bad  plate  data;  check  data  deck. 

(6)  NEGATIVE  MOMENT  OF  INERTIA,  PLATE  NO.  XX— Indicates  negative 
moment  of  inertia  for  plate;  program  converts  value  to  positive. 

(7)  NEGATIVE  THICKNESS,  PLATE  NO.  XX— Indicates  negative  tQ,  t^,  ty, 

or  t  ;  program  converts  value  to  positive, 
s 

(8)  QUADRILATERAL  PLATE  NO.  XX — Has  a  re-entrant  corner.  Local 
coordinates  are:  X  coordinates  (XX,  XX,  XX);  Y  coordinates 
(XX,  XX,  XX).  Indicates  plate  has  an  interior  angle  greater  than  180 
degrees;  if  data  are  valid,  then  represent  the  quadrilateral  as  two  triangles. 
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(a)  Positive  In-Plane  Plate  Stresses  Per  Linear  Inch  At  CG  Of  Plate 


(b)  Positive  Out-Of-Plane  Plate  Stresses  Per  Linear  Inch  At  CG  Of  Plate 


Figure  13.  Plate  Stresses 
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(a)  Positive  Beam  Stresses-Straight  Beams 

Y 


Figure  14.  Beam  Stresses 
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(9)  QUADRILATERAL  PLATE  NO.  XX/DEGREE  OF  NONCOPLA¬ 
NARITY  XX  PERCENT/ACCEPTABLE-ANALYSIS  CONTINUED 
(larger  than  0. 1  percent)  —  Indicates  the  four  plate  nodes  are 
not  coplanor,  but  by  an  acceptable  amount;  check  coordinates 
of  the  nodes. 

(10)  QUADRILATE RIAL  PLATE  NO.  XX/DEGREE  OF  NONCOPLA¬ 
NARITY  XX  PERCENT/NOT  ACCEPTABLE -ANALYSIS  DELETED— 
Indicates  that  the  four  plate  nodes  are  not  coplanor  by  an 
unsatisfactory  amount;  check  coordinates  of  the  nodes. 

(11)  INCORRECT  NODE  NUMBER,  BEAM  NO.  XX  —  Indicates 
beam's  node  number  is  greater  than  the  number  of  nodes  in  the 
structure. 

(12)  NEGATIVE  MODULUS  OF  ELASTICITY,  BEAM  NO.  XX  — 

Indicates  negative  Young's  Modulus  for  this  beam;  program 
converts  value  to  positive. 

(13)  NEGATIVE  MODULUS  OF  RIGIDITY,  BEAM  NO.  XX  —  Indicates 
negative  shear  modulus  for  this  beam;  program  converts  value 

to  positive. 

(14)  INPUT  NOT  IN  PROPER  SEQUENCE,  BEAM  NO.  XX  SHOULD  BE 
XX  —  Indicates  bad  beam  data;  check  data  deck. 

(15)  SECTION  PROPERTY  XX  IS  OUT  OF  SEQUENCE  —  Indicates 
section  properties  for  this  beam  are  not  sequential;  re-order 
section  property  data. 

(16)  THE  NUMBER  OF  PARTITIONS  IN  K(FF)  HAS  BECOME  GREATER 
THAN  800  IN  THE  UNKNOWN  DEFLECTION  COMPUTATIONS  — 
Indicates  too  many  stiffness  partitions  have  been  created;  reduce 
size  of  structure,  if  possible. 

(17)  MATRIX  IS  SINGULAR  —  Indicates  a  diagonal  K  partition  is 
singular;  check  data  (Either  the  structure  is  unstable,  or  some 
orientation  nodes  are  not  constrained. ) 

The  following  diagnostic  messages  from  the  eigenvalue  module  indicate 

errors  in  reading  and  writing  tape. 

(1)  FLEXIBILITY  MATRIX  COULD  NOT  BE  FOUND 
CURRENT  TAPE  IN  USE  =  XX,  ERROR  CODE  =  XX 


(2)  SPACING  ERROR  OCCURRED  WHILE  TRYING  TO  WRITE  ON  BINARY  TAPE. 
NFILE  =  XX,  NMAT  =  XX 

CURRENT  TAPE  IN  USE  =  XX.  ERROR  CODE  =  XX 

(3)  ERROR  WHILE  j  READING  MASS  MATRIX  l  ON  TAPE 

(WRITING  DYNAMIC  MATRIX  f 

CURRENT  TAPE  IN  USE  =  XX,  ERROR  CODE  =  XX. 

In  the  preceding  diagnostics,  the  logical  tape  number  is  given,  NFILE  is  the 
number  of  file  marks,  and  NMAT  is  the  number  of  matrices.  Error  codes 
are  as  follows: 

(1)  File  spacing  error 

(2)  Negative  matrix  spacing 

(3)  Matrix  spacing  error 

(4)  Check  sum  error 

(5)  Tape  name  wrong 

Other  possible  error  messages  are: 

(1)  EIGENVALUE  =  XX  FOR  MODE  XX,  HAS  A  NEGATIVE  VALUE,  WHICH 
DOES  NOT  DESCRIBE  A  PROPER  PHYSICAL  SYSTEM. 

(2)  MODE  HAD  ONLY  XX  PERCENT  ACCURACY,  WHICH  IS  UNDER  THE 
GIVEN  LIMIT  OF  XX.  THE  TWO  VECTORS  WHICH  SHOULD  BE 
IDENTICAL  ARE  GIVEN  BELOW. 

XXXXXX  1  xxxxxx 

2 

3 
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SAMPLE  PROBLEMS 

This  section  describes  the  application  of  RANVIB  to  random  structural 
response  analysis.  Example  problems  are  given  for  two  structures:  a  simply 
supported  beam  and  a  sandwich  plate.  System  input  and  output  data  described 
in  this  section  are  given  in  appendix  III. 

1.  SIMPLY  SUPPORTED  BEAM 

The  first  example  illustrates  solution  options  of  the  program  on  a  very 
simple  structure.  A  simply  supported  beam  is  shown  schematically  in  figure  15 
together  with  its  properties  and  structural  idealization.  The  beam  is  loaded  by 
a  random  pressure  wave  convected  along  the  beam  in  the  x  direction;  its  charac¬ 
teristics  are  shown  in  the  figure.  Deflection  and  stress  responses  of  the  beam 
are  sought  for  this  loading. 

Input  data  for  phase  I  are  given  in  appendix  III,  page  85.  Columns  79  and 
80  are  used  in  these  examples  to  identify  sections  and  cards  defined  in  the 
input-preparation  description.  The  beam  is  idealized  into  six  equal  elements, 
and  beam  bending  is  the  only  deformation  considered.  Pressure  is  applied 
normal  to  the  beam’s  surface,  and  rotary  inertia  effects  are  assumed  to  be  neg¬ 
ligible.  Therefore,  only  transverse  deflections  are  retained  explicitly  in  the 
analysis.  Four  mode  shapes  and  frequencies  are  to  be  determined  in  this  example. 

Printed  results  of  the  phase  I  analysis  are  presented  in  appendix  III,  pages 
86  through  96. 

Several  phase  II  solution  options  are  presented  for  the  beam  problem 
to  illustrate  their  use.  The  analyses  performed  are:  (1)  option  3  joint  moments, 
(2)  option  2  joint  moments,  (3)  option  3  cross  PSD,  (4)  option  2  cross  PSD, 
and  (5)  option  1  cross  PSD. 

Input  data  for  the  option  3  joint -moment  solution  are  given  in  appendix 
III,  page  97.  Deflection  and  stress  covariance  and  second  spectral  moment 
matrices  are  to  be  computed.  Structural  damping  is  assumed  with  a  damping 
coefficient  of  0.  010.  Thus,  the  critical  camping  ratio  is  0.  005  for  all  freq¬ 
uencies  considered.  The  force  cross -PSD  matrix  is  printed  only  for  the  first 
natuial  frequency.  Coordinates  across  the  width  of  the  beam  were  selected 
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ct 


1  2  3  4 _ 5  6  7 


6x9  =  54  IN. 


PRESSURE  WAVE  VELOCITY . ct  =  7,200  IPS 

PRESSURE-WAVE  DECAY  CONSTANT . d  =  0.001 189  IN.'1 

YOUNG'S  MODULUS . E  =  10.5  x  10b  PS  I 

MOMENT  OF  INERTIA  I . I  =  0.00327  IN.4 

WEIGHT . V  =  0.101  LB/IN. 


Figure  15.  Simply  Supported  Beam 

such  that  the  nodal  areas  calculated  are  in  unity.  Pressure  PSD  is  assumed 
to  be  in  unity  at  each  of  the  four  natural  frequencies  used  in  the  analysis. 

Analysis  results  are  presented  in  appendix  III,  page  98.  Diagonal  terms 
of  the  deflection  covariance  matrix  are  mean-square  values  of  the  modal 
deflections.  In  the  stress  covariance  matrices,  values  are  obtained  only 
for  Mv  and  F„  components,  because  twisting,  axial,  and  out-of-plane  de for- 
mations  are  not  considered.  Elements  of  the  second  spectral  moment  matrices 
can  be  used  in  equation  (47)  to  predict  the  number  of  zero  crossings  of  response 
time  history.  As  an  example,  this  number  for  the  shearing  force  F  at  the 

end  of  the  beam  is 

„  „  1  r  2.316  x  10°  1  1/2 

°  17  L  3.^98  x  101*  J 

=  25.9/seC 

Thus,  the  "apparent"  frequency  is 

\  N  =  13  Hz 

2  o 
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For  the  moment  at  the  beam  center,  this  frequency  is 


1 

2 


N 

o 


1  [~  1.801  X  1Q10~|  -1-/2 

2ir  L  1.046  X  10T  J 


6.6  Hz 


Input  data  for  the  option  2  joint-moments  analysis  are  shown  in  appendix  in, 


page  111.  One  cross  product  for  each  mode  is  to  be  included  in  the  calculations. 

/  GJ-i  +  CJo  Clin  +  0t)o 

Pressure  PSD  for  each  of  the  cross-term  frequencies  ( — — — ,  — ,  and 

Wo  +  UJ4  \ 

■  ---0 - )  is  also  assumed  to  be  unity. 


Printed  analysis  results  are  presented  in  appendix  IE,  page  112.  In  add¬ 
ition  to  force  cross  PSD  being  printed  for  the  first  modal  frequency,  force 
co-PSD  and  quad-PSD  matrices  are  printed  for  the  first  cross  modal  term.  The 
real  part  of  the  response  is  not  significantly  different  from  that  obtained  in  option 
3.  For  this  problem  in  which  damping  is  very  low  and  modal  frequencies  are 
weU  separated,  the  cross  terms  do  not  contribute  significantly  to  the  result. 
Imaginary  parts  of  the  joint -moment  matrices  indicate  phase  differences 
between  responses. 

An  illustration  of  input  data  for  option  3  cross -PSD  calculation  is  given  in 
appendix  in,  page  137.  Both  deflection  and  stress  cross  PSD’s  are  to  be 
computed  for  one  frequency:  the  first  natural  frequency  of  the  beam.  Output 
results  are  presented  in  appendix  EH,  page  138. 

Input  data  for  an  option  2  analysis  for  the  preceding  problem  is  given  in 
appendix  IH,  page  148.  Printed  results  for  this  analysis  are  presented  in 
appendix  m,  page  149. 

The  final  example  given  for  the  beam  problem  is  an  option  1  cross -PSD 
analysis  for  the  same  frequency  used  in  the  previous  solutions .  Input  data  are 
shown  in  appendix  IH,  page  166.  For  this  illustrative  problem,  the  matrix 
of  viscous  damping  coefficients  has  been  selected  to  be  proportional  to  the 
stiffness  matrix.  Normally,  option  1  would  not  be  used  in  this  case,  but 
these  coefficients  were  chosen  so  that  the  damping  ratio  at  the  first  natural  fre¬ 
quency  is  equal  to  that  used  in  previous  solutions.  Thus,  the  results  of  this 
analysis  can  be  compared  to  the  other  options . 
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Printed  results  of  the  analysis  are  presented  in  appendix  III,  page  167. 
Computed  response  is  essentially  the  same  as  in  other  options.  If  joint 
statistical  moments  are  computed  in  option  1,  the  response  spectrum  must 
be  determined  at  a  sufficient  number  of  frequencies  since  a  numerical  integra¬ 
tion  is  performed  to  calculate  joint  moments . 

2.  SANDWICH  PLATE 

The  second  example  problem  is  for  a  sandwich  plate  simply  supported  on 
two  opposite  sides  and  supported  on  beams  on  the  other  two  sides.  This 
example  is  offered  as  an  illustration  of  phase  I  input  data  for  more-complex 
structures. 

The  structure  is  shown  schematically  in  figure  16.  Plate  and  beam 
properties  are  also  shown  in  the  figure.  The  neutral  axis  of  each  beam  is  1 
inch  below  the  plane  of  the  plate. 

Structural  idealization  used  in  the  problem  is  shown  in  figure  17.  Input 
data  for  phase  I  are  presented  in  appendix  III,  page  183.  All  in-plane  deform¬ 
ations  are  assumed  to  be  negligible  and  ,  thus,  are  not  considered.  Shear 
deformation  of  the  sandwich  core  material  is  also  neglected.  Only  bending  and 
torsional  distortions  are  considered  for  the  beams.  Transverse  deflections  at 
the  nodes  not  rigidly  supported  are  retained  in  the  analysis.  Four  modal  fre¬ 
quencies  are  to  be  determined. 

The  printed  phase  I  output  for  the  problem  is  presented  in  appendix  HI, 
page  186. 


68 


24  IN. 


PLATE  PROPERTIES: 


YOUNG'S  MODULUS  .  . 
POISSON'S  RATIO  .  .  . 
PLATE  (SKIN)  THICKNESS 
SANDWICH  DEPTH  .  . 
WEIGHT . 


E  =  10.5  x  106  PSI 
v  =  0.33 
lo/2  =  0.02  IN. 
d  =  0.5  IN. 
y  =  o.oo7  psi 


BEAM  PROPERTIES: 


YOUNG'S  MODULUS  . 
SHEAR  MODULUS .  . 
MOMENT  OF  INERTIA 
TORSION  CONSTANT 
WEIGHT . 


E  =  10.5  x  IQ6  PSI 
G  =  3.9  x  106  PSI 
IY  =  0.1  IN.4 
J  =  8  x  10‘5  IN.4 
y  =  0.016  LB/IN. 


Figure  16.  Sandwich  Plate 


69 


y 


NODE  NUMBER 
PLATE  NUMBER 

BEAM  NUMBER 


Figure  1 7.  Idealization  of  Sandwich  Plate 
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APPENDIX  I 


STRUCTURAL  ELEMENTS 

1.  A  BRIEF  DISCUSSION  OF  ELEMENT  STIFFNESS  MATRICES 

As  implied  in  the  input  description,  the  structural  matrix  module  has  all 
of  the  standard  elements  of  any  displacement  program,  i.e. ,  the  rod,  the  beam, 
and  a  plate  capable  of  carrying  in-plane  and/or  out-of-plane  loading. 

The  element  flexibility  matrix  for  the  curved  beam  is  obtained  using  the 
principle  of  complementary  energy.  This  matrix  is  then  inverted  to  give  the 
element  stiffness  matrix.  The  element  stiffness  matrix  for  the  straight  beam 
is  obtained  from  that  for  the  curved  by  considering  the  limit  as  the  radius 
approaches  infinity.  The  rod  stiffness  is  simply  that  for  the  straight  beam, 
assuming  no  bending  resistance. 

Matrices  for  quadrilateral  plate  elements  are  generated  by  first  subdividing 
the  elements  into  four  triangles . 

The  triangular  plate  element  is  basically  that  found  in  reference  9.  Out- 
of-plane  bending  stiffness  is  obtained  by  "sandwiching"  two  identical  in-plane 
plates  together.  The  plate  does  have  the  unique  characteristic  of  having  an 
in-plane  moment  stiffness.  This  is  used  to  prevent  singularities  when  a  user 
does  not  wish  to  support  this  freedom. 

For  a  detailed  description  of  the  elements  in  this  program,  see  reference  4. 

2.  GEOMETRIC  STIFFNESS  MATRICES 

The  program  contains  geometric  stiffness  matrices  for  the  following: 

a.  Beam  —  axial  load  in  local  X  direction  for  rotation  about  the  local  Y  axis, 
or  rotation  about  the  local  Z  axis 

b.  Triangular  plate  —  in -plane  stress  in  the  local  X  direction 

c.  Triangular  plate  —  in-plane  stress  in  the  local  Y  direction 

d.  Triangular  plate  —  in-plane  pure  shear  stress 

For  a  description  of  the  use  and  derivation  of  these  elemental  matrices, 
see  reference  6. 
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Figure  18  shows  the  matrix  elements  for  a  beam  of  length  L  with  tensile  force 
Fx  and  rotation  about  the  local  Y  axis.  This  matrix  can  be  modified  for  rota¬ 
tion  about  the  local  Z  axis  by  shifting  rows  and  columns  3,  5,  9,  and  11  to  rows 
and  columns,  2,  6,  8,  and  12,  respectively. 

Geometric  stiffness  matrices  for  a  triangular  plate,  such  as  those  shown  in 
figure  7(b),  with  in-plane  tensile  stresses  a  and  a  and  shear  stress  t  are 

-A.  I  A  Y 

shown  in  figure  19.  The  thickness  of  the  plate  is  t  and  XQO  =  X  -  X 

32  3  2 

3.  INTERPRETATION  OF  PLATE  INPUT  DATA 

For  a  uniform  and  solid  plate,  the  actual  thickness  t  is  read  from  the 
plate  card.  Other  plate  properties  are  computed  by  the  program  as  follows: 

tx  -  °  *x  =  0 

*Y  =  0  4s  =  4o 

tv  =  0  I  =  t  3/l2 

I  so 

The  values  of  tQ  and  t  must  be  positive.  For  negative  values,  the  pro¬ 
gram  takes  the  absolute  value  and  prints  a  warning. 

For  stiffened  or  sandwich  plates,  values  for  the  properties  in  addition  to 

t  are  read  from  card  3  of  the  plate  data .  Stiffener  properties  t  ,  I  .  ,  t 
u  X  Y  Y  ’ 

and  L.  can  be  zero.  If  the  effective  shear  thickness  t  is  read  in  as  zero,  the 
-X  s 

program  sets  t  equal  to  t  .If  the  plate  bending  moment  of  inertia  I  is 

S  O  g  S 

read  in  as  zero,  the  program  sets  I  =  t  /12  .  The  quantities  t  ,  t 

so  o  X 

ly  »  *x  »  *s  ’  and  Js  must  not  be  negative.  For  negative  values ,  the  pro¬ 
gram  takes  the  absolute  value  and  prints  a  warning. 

The  stiffened-plate  option  is  one  alternative  when  structure  can  be 
idealized  into  an  assembly  of  equivalent  orthotropic  plate  elements.  The  model 
for  this  plate  element  uniformly  distributes  equivalent  stiffener  properties 
into  the  plate.  An  alternate  procedure  with  better  results  when  stiffeners 
are  widely  separated  is  to  use  beam  elements  to  represent  the  stiffeners. 
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Figure  18.  Geometric  Stiffness  Matrix  for  a  Beam 
(Rotation  about  Local  Y  Axis) 
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Figure  19.  Geometric  Stiffness  Matrices  for  a  Triangular  Plate 
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APPENDIX  H 


ADMITTANCE  INTEGRALS 


To  obtain  joint  response  moments  in  the  normal  mode  analysis  (option  2), 
it  is  necessary  to  evaluate  certain  integrals  having  system  admittance  terms. 
These  integrals  have  been  evaluated  in  reference  8,  and  results  are  presented 
in  this  appendix. 

The  admittance  integrals  appearing  in  the  equation  for  deflection 
covariance,  equation  (63),  page  27 ,  are  of  the  form: 


0 


(D.D.  +  E.E, )  dw 
1  j  x  j 


oo 

0 


D.  E.  dw 
i  J 


where  D-  and  E.  are  defined  by  equations  (60)  and  (61),  page  26.  For 
J  J 

convenience  in  expressing  the  integrals,  the  modal  damping  is  specified 
by  a  parameter  \i  j  given  as 

y.  =  2 e.w. 

j  j  j 

Damping  is  assumed  to  be  small  (fj<  1)  •  Results  of  the  integrations  are 

oo 

(  D.D.  +  E.E. )  dw  = 

'  x  j  i  J 


tt(  +  y  j ) 


NLM.  Luk**  +  (y j2  +  “i2  +  “j4  +  +  2u)± 


,2] 
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rOO 


D.  E.  dw  = 
i  J 


—  (l—  loSe  — 4  +  Uij  -  —  (^i2  '  2a)  2)]/_ 

2M.M.  I  2  w.4  J  2 

i  J  i  ill 


<Au>.2  -  M.2 


B. 


+  [c  .  +  —  (y.2  -  2m  2)] 

J  ^  J  J 


(  M  .  \/  4o^2  -  Mj2  )  | 


where  the  constants  are  given  by 


(o)  4-  0)  4)  [-y.o)2  (y,Z-2o)  2)-  m.0).4]+m.0).2[0).  4  (m.2-2oj,2)-  w.4  (m.2-2w.2)] 
J  J-*--1--*-  J-*-  J  J  J  J  1  1 

J.4[m.  2-y  .2-2(to  2-w  2  )  ]  [o) .  4  (  m  •  2-2w  2  )  -  ui .  4(  m  .  2— 2d).2  )  ]  -  u)  4  (w.  4-o)  4)2 

«J  J  J  J  1  1  1  I  J 


B.  .  = 
ij 


A.  . 
1J 


Bii  “i  ~  y,2  ~  2  (ui  2  -  w.2)]  -  m  .0). 2 

-1-  J  1  J  J  1 


0).  4  -  O)  .  4 

1  J 


y  .  w.  2  -  ui .  4  A.  . 
J  1  J  iJ 


C.  .  = 

iJ 


0) . 
l 


The  expression  for  the  second  integral  is  valid  only  for  i  =  j  .  This 
integral  is  not  required  for  i  =  j  ,  because  it  is  multiplied  by  a  null  matrix 
in  equation  (63),  page  27. 

Calculation  of  second  spectral  moments  requires  the  evaluation  of  integrals 

2 

formed  from  the  previous  ones  by  multiplying  the  integrands  by  cm  .  These 
integrals  are  given  by 


u)2  (D.  D.  +  E.  E.)  do)  — 
1  J  1  J 


TT  (  M  -  0).2  +  y.  o).2) 
1  J  J  1 


[M.M.]  [<jj.  4  +  (m.2  +  y .  M . )  o).2  +  O).4  +  (m-2  +  y-y.)  O).2  -  2w.2u).2] 
1  j  1  j  1  J  1  J  1  1  J  J  1  J 
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CO 


D .  E .  dxo  = 
i  J 


,  (  F.  .  oo .  F.  .  u 

— i -  |  1J  lop  _J — +  fE. - (u.* 2  -  2o). 2 )  1  - 

2M-  M-  I  2  6  4  U  2  ^x  X  \Tkii}.  2-y .  2 


1  J 


U)  . 
1 


1  11 


F.  . 


+  [G.  .  +  (y  .2  -  2oo  2)] 
xj  2  j  J 


y.  v Uo) . 2-y  . 2 
J  J  J 


where 


y  a)  6(o>  **-u)  4)  +  y  w  4[<o  4  ( y .  2-2u)  . 2 )  -  a,.4  (y.2-2(D.  2)] 

«-)  ^ _ J _ J  J _ 1 _ J _ J _ 1 _ 1 _ 1 

.  4  [  y ,  2- y ^  2-2 ( (D.  2-^  2  )  ]  [  0).  ^  2-2(1) 2  ) -0) .  4  (  y .  2 -2a, .  2  )  ] -d).  4  ( (D^^  ^-O) .  ^  ) 2 


F.  .  = 

1J  oo.  4  Ti.  2—i.  2 


E.  .  = 

xj 


F.  .  a).4  [y  .  z-y  .  2-2(  a,.  2-0) .  2  )  ]  -  y.u>.4 
x  J  l  l  J  i  J_ J  x 

4  4 

OJ.  ~  03  . 


G.  .  =  -  E.  .  u 

XJ  U  4 

l 


Again,  it's  not  necessary  to  evaluate  the  second  of  these  integrals  for 


i  =  J  • 
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SYSTEM  INPUT  AND  OUTPUT  DATA 
The  following  is  a  list  of  the  information  in  this  appendix: 

Page 

Phase  I  System  Control  Cards .  80 

Phase  II  System  Control  Cards .  83 

Phase  I  Input  Data — Beam .  85 

Phase  I  Output  Data — Beam  .  86 

Phase  II  Input  Data — Option  3,  Joint  Moments  .  97 

Phase  II  Output  Data — Option  3,  Joint  Moments .  98 

Phase  II  Input  Data — Option  2,  Joint  Moments  . Ill 

Phase  II  Output  Data — Option  2,  Joint  Moments . 112 

Phase  II  Input  Data — Option  3,  Cross  PSD . 137 

Phase  II  Output  Data — Option  3,  Cross  PSD . 138 

Phase  II  Input  Data — Option  2,  Cross  PSD . 148 

Phase  II  Output  Data — Option  2,  Cross  PSD  . 149 

Phase  II  Input  Data — Option  1,  Cross  PSD . 186 

Phase  II  Output  Data — Option  1,  Cross  PSD . . 167 

Phase  I  Input  Data — Sandwich  Plate . 183 

Phase  I  Output  Data — Sandwich  Plate . 186 
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PHASE  1  SYSTEM  CONTROL  CARDS 


SEXECUTE 

IBJOB 

$ IBJOB 

SPOOL 

BUFCT=07* 'UNIT02' # 

SETC 

•  UN I T 1 1 • » 'UNIT12' i 

SETC 

•  UNI T 17 • *  »UNIT10« » 

SGROUP 

OPNCT=07»BUFCT=07» 

SETC 

•UNITll • »  'UNIT12*  » 

SETC 

•  UN  I T 1 7 • ♦  'UNIT10' » 

SIEDIT 

SYSLB4 » SRCH 

SIBLDR 

BLK 

SIBLDR 

PHASE  1 

SIBLDR 

WRTETP 

SIBLDR 

READTP 

SIBLDR 

FVIO 

SIBLDR 

FILE 

SIBLDR 

FBSF 

SIBLDR 

FFSF 

SIBLDR 

FFSR 

SIBLDR 

fkrd 

SIBLDR 

FRUN 

SORIGIN 

PHASE1 »SYSCK2 

SIBLDR 

MAST* 

SIBLDR 

PAGH* 

SIBLDR 

UNPAC* 

SIBLDR 

PRINT* 

SORIGIN 

ALPHA  1 » SYSCK2 

SIBLDR 

SUBM1* 

SORIGIN 

BETA1 »SYSCK2 

SIBLDR 

GENRA* 

SIBLDR 

REDUC* 

SORIGIN 

GAMMA  1 »SYSCK2 

SIBLDR 

INFO* 

SORIGIN 

GAMMA1 *  SYSCK2 

SIBLDR 

PLATE* 

SIBLDR 

MUL  1* 

SIBLDR 

MUL2* 

SIBLDR 

PST  IF* 

SIBLDR 

QUAD* 

SIBLDR 

LAMK* 

SIBLDR 

KLAMT* 

SIBLDR 

TRI* 

$  I BLDR 

INP* 

SIBLDR 

INPM* 

SIBLDR 

INPST* 

SIBLDR 

OUT  P* 

SIBLDR 

OUT  PM* 

SIBLDR 

OUTPS* 

SIBLDR 

COMBI* 

SIBLDR 

STORE* 

SIBLDR 

MOVE* 

SIBLDR 

PMTR* 

•UN 
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•  UN 
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•  UN 
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* • UN  IT  16* 

•UN 

I  T03  ’ 
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$  I BLDR  SMTR* 

$  I BLDR  LOCAL* 
$  I BLDR  COPLA* 
SORIGIN 
S INCLUDE 
$  I BLDR  BEAM* 

$  I BLDR  TINVR* 
$  I BLDR  SMULT* 
$  I BLDR  MULT* 

5  I BLDR  SBMTR* 
$  I BLDR  SSTIF* 
$  I BLDR  MAD* 

$  I BLDR  SBGS* 

S I BLDR  OFST * 

$  I BLDR  CSTIF* 

$  I BLDR  CBMTR* 

SORIGIN 

$  I BLDR  MERGE* 

$  I BLDR  SOR* 

SORIGIN 

$  I BLDR  MERGE* 

SORIGIN 

S I BLDR  STRES* 

S I BLDR  MSTRE* 

SORIGIN 

S I BLDR  SORCON 

$  I BLDR  TEST* 

SORIGIN 

$  I  BLDR  FKSRT 

SORIGIN 

$  I BLDR  KFSRT 

SORIGIN 

$  I BLDR  CNCT* 

SORIGIN 

S I BLDR  EXPND* 

SORIGIN 

$  I BLDR  EXTRN* 

SORIGIN 

$  I  BLDR  DELET* 

S  I  BLDR  SSORT* 

SORIGIN 

S  I  BLDR  HELP 

$  I  BLDR  YTLOSB 

S  I  BLDR  INVERT 

$  I BLDR  DAT ASB 

SORIGIN 

$  I BLDR  FREMO* 

$  I BLDR  WDVALV 

S I BLDR  INRPRD 

SORIGIN 

$  I BLDR  WDHESS 

SORIGIN 


GAMMA  1 »SYSCK2 
FSCN#  FXP2 


BETA1 #SYSCK2 

DELTA1 »  SYSCK2 
DELTA1 »  SYSCK2 

ALPHA  1 » SYSCK2 

SORT1 1  SYSCK.2 
SORTl »  SYSCK2 
SORT1 #SYSCK2 
EXP  1  *  SYSCK2 
EXP  1 »SYSCK2 
SORTl #SYSCK2 

ALPHA 1 * SYSCK2 


PHASE 1 #SYSCK2 

ALPHA  »SYSCK2 
ALPHA » S YSCK2 
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$  I BLDR  WDQRIT 

SORIGIN 

$  I BLDR  WDSORT 

SORIGIN 

$  I BLDR  WDVECT 

SORIGIN 

S  I  BLDR  WDTRN1 

SORIGIN 

$  I BLDR  WDTRN2 

SORIGIN 

$  I BLDR  AMERG* 

SORIGIN 

S I BLDR  SMERG* 


ALPHA  »SYSCK2 
ALPHA  »SYSCK2 
ALPHA  »S YSCK2 
ALPHA  »SYSCK2 
PHASE1 *SYSCK2 
PHASE1 *  SYSCK2 
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PHASE 

II  SYSTEM  CONTROL  CARDS 

SEXECUTE 

IB  JOB 

SIBJOB 

SPOOL 

BUFCT=09* 'UNIT01*  * 

SETC 

•UNI T08 • •  MJNIT10* • 

SETC 

•  UNIT  15 1 # IUNIT16*  * 

SGROUP 

OPNCT=09* 'UNITOI* * 

SETC 

•UNIT08 • » 'UNITIO' » 

SETC 

•  UNI T 1 5 • *  «UNIT16« • 

SIEDIT 

SYSLB4.SRCH 

SIBLDR 

PHASE2 

SIBLDR 

SCALE* 

SIBLDR 

FVI 0 

SIBLDR 

FILE 

SIBLDR 

READTP 

SIBLDR 

WRTETP 

SIBLDR 

FBSF 

SIBLDR 

FFSF 

SIBLDR 

FFSR 

SIBLDR 

FKRD 

SIBLDR 

FRUN 

SORIGIN 

ALPHA  t  SYSCK.2 

SIBLDR 

HELP 

SIBLDR 

YTLOSB 

SIBLDR 

INVERT 

SIBLDR 

DATASB 

SORIGIN 

ALPHA  #  SYSCK2 

SIBLDR 

TAPOS* 

SORIGIN 

ALPHA »SYSCK2 

SIBLDR 

RANLO* 

SORIGIN 

BETA.SYSCK2 

SIBLDR 

ARIA** 

SORIGIN 

BETA.SYSCK2 

SIBLDR 

CONST* 

SORIGIN 

BETA  *SYSCK2 

SIBLDR 

NO I SO* 

SIBLDR 

OUTPUT 

SORIGIN 

ALPHA  »SYSCK2 

SIBLDR 

PEDAN* 
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ALPHA#SYSCK2 

SIBLDR 

PRNTA* 

SORIGIN 

ALPHA  #  SYSCK2 

SIBLDR 

PRNTB* 

SORIGIN 

ALPHA  #SYSCK2 

SIBLDR 

CONS* 

SORIGIN 

ALPHA *SYSCK2 

SIBLDR 

PRNTC* 

UN 
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» 

UN 
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UN 
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UN 
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» • UN  I T04 • 

t 

UN 
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» 'UNIT  12' i 

, 1  UN  I T 13 ' 

UNI T 14 • 

UN 

I  T17  • 
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SORIGIN 

SIBLDR  PRNTD* 

SORIGIN 

$  I BLDR  DSEC1* 

SORIGIN 

S I BLDR  ADMN3* 

SORIGIN 

$  I BLDR  ADDMA* 

SORIGIN 

$  I BLDR  ADMN2* 

SORIGIN 

$  I BLDR  CQJD* 

SORIGIN 

$  I BLDR  SUMT* 

SORIGIN 

$  I BLDR  SUM2* 

SORIGIN 

$  I BLDR  SECM3* 

SORIGIN 

$  I BLDR  SECM2* 

SORIGIN 

$  I BLDR  ADMT3* 

SORIGIN 

$  I BLDR  ADMT2* 

SORIGIN 

$  I  BLDR  CQCPS* 

SORIGIN 

S I BLDR  SUM3* 

SORIGIN 

$  I BLDR  PRNTE* 

SORIGIN 

$1 BLDR  DSECM* 


ALPHA  »SYSCK2 
ALPHA  »SYSCK2 
ALPHA  »  SYSCK2 
ALPHA  #  SYSCK.2 
ALPHA  »SYSCK2 
ALPHA ♦  SYSCK.  2 

ALPHA » SYSCK2 
ALPHA  #SYSCK2 
ALPHA i SYSCK2 
ALPHA  *SYSCK2 
ALPHA  #SYSCK.2 
ALPHA  » SYSCK2 
ALPHA  »SYSCK2 
ALPHA  »SYSCK2 
ALPHA « S YSCK2 
ALPHA ♦ SYSCK2 


84 


PHASE  I  INPUT  DATA  -  BEAM 
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PHASE  I  OUTPUT  DATA  --  BEAM 
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SIMPLY  SUPPORTED  BEAM 
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SIMPLY  SUPPORTED  BEAM 
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STIFFNESS  MATRIX 
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BELOW  IS  THE  COMPLETE  FLEXIBILITY  MATRIX  UNSCALED.  THE  SCALE  FACTOR 
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BELOW  IS  THE  DYNAMIC  MATRIX 
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BELOW  IS  THE  INERTIA  MATRIX 


* 


o  o  o  o 

-O  CD  "O  00 

in  o  m  o 

O  O  sj- 

•  O'  ^  o 

o  o  o  o 
o  o  o  o 

•  •  •  • 

o  o  o  o 


<  IN 


>}• 


94 


BELOW  ARE  TEE  SHAPES  FOR  MODES 
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TABLE  CF  FINAL  RESULTS 
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PHASE  II  INPUT  DATA -OPTION  3,  JO  I  NT  MOMENTS 
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PHASE  II  OUTPUT  DATA  --  OPTION  3,  JOINT  IV\OIV\ENTS 
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The  following  force  cross  power  spectral  oensity  matrices  were  generated 
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PHASE  1 1  OUTPUT  DATA  -  OPTION  2,  CROSS  PSD 
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FORCE  CROSS  POWER  SPECTRAL  DENSITY  FOR  DECAYED  PROGRESSIVE  WAVES  FOR  USE  WITH  PANEL 
SIMPLY  SUPPORTED  BEAM,  OPTION  2  CROSS  PSD  B1 


3 

a 

Z 

p- 

CM 

o 

O 

or 

< 

• 

UJ 

U 

O 

z 

o 

Ui 

o 

o 

a: 

•— * 

o 

UJ 

p- 

o 

X 

u 

m 

UJ 

* 

>- 

CL 

• 

UI 

o 

X 

CD 

P- 

O 

co 

X 

z 

CM 

o 

m 

u 

M 

O 

o 

1 

OC 

UJ 

X 

UJ 

UJ 

UJ 

o 

X 

co 

O 

o 

o 

X 

o 

o* 

o 

z 

3 

o 

0* 

o 

< 

z 

1 

o 

o 

o 

o 

O' 

o 

p- 

fit 

UJ 

UJ 

o 

O' 

o 

UJ 

o 

oc 

O' 

o 

< 

o 

m  o 

< 

• 

• 

• 

- 

CL 

z 

o 

o 

o 

o 

a 

z 

o 

o 

■ 

UJ 

1  -• 

►u 

p- 

X 

-J 

p- 

CM 

o 

•< 

P- 

o  g 

o 

o 

o 

p- 

3 

z  v 

UJ 

UJ 

z 

—•  u 

CL 

UJ 

UJ 

X 

d 

OC 

o 

o 

p- 

UJ  UJ 

1 

CD 

o 

O' 

Z 

CL 

co  X 

o 

O' 

» 

< 

X  P- 

UJ 

g 

o 

o* 

CO 

UJ 

3 

CL 

•— « 

o 

O' 

M 

Q. 

z  o 

< 

O' 

X 

a. 

p- 

g 

CM 

O' 

< 

< 

UJ  • 

z 

>* 

• 

• 

♦  1 

O  OC 

o 

g 

O 

o 

O  X 

>- 

o  < 

*-• 

UJ 

Z  _J 

UJ 

UJ  Q 

X 

3 

g 

X 

o  z 

p- 

0-t 

o  o 

UJ 

o 

p- 

CM 

o 

o  < 

(-4 

oc 

o 

O 

o 

o 

-J  Q 

M 

o 

o  z 

Q 

■ 

o  z 

o 

UJ 

p- 

UJ 

UJ 

o  o 

UJ 

>■  UJ 

o 

o 

o 

o  — 

P- 

CM  CD 

O  a. 

o 

CM 

OC 

o 

O' 

CM  P- 

U 

< 

oon 

oc 

CM 

< 

o 

O' 

r-  o 

UJ 

w 

1  1 

U.  UJ 

-J 

CM 

_J 

O 

1 

O' 

•  UJ 

-J 

UJ  UJ  UJ 

O  a. 

o 

CM 

3 

O 

o* 

O  ac 

UJ 

o  o  o 

> 

CM 

o 

CO 

UJ 

O' 

• 

«/> 

O 

o  o  o 

z  z 

g 

<M 

•— < 

•-» 

Q 

O' 

o 

o  o  o 

o  o 

• 

o 

• 

O 

• 

UI 

■ 

z 

o  o  o 

UJ 

O 

z 

O 

z 

o 

ac 

z 

UJ 

■ 

<7*  O  O 

p"  p- 

X 

UJ 

< 

H  O 

UJ 

CO  o  o 

g  o 

P- 

CL 

o 

O  M 

CO 

-M  CM  o 

UJ  UJ 

OC 

UI 

«/) 

p- 

n 

r*-  o 

OC  oc 

z 

o 

UI 

Z 

o 

UJ 

•  < 

UJ 

w 

•  •  • 

•-*  M 

o 

o. 

o 

o 

UI  <3 

>■ 

o  o  o 

O  CD 

«* 

p- 

>  < 

< 

1 

CO 

UJ 

«/> 

UJ 

P- 

UJ 

•— 

<  Cl 

X 

cm 

UJ  UJ 

UJ 

O 

UJ 

o 

UJ 

o 

Ui 

• 

X  o 

■  ■  « 

X  X 

g 

—4 

o 

o 

ac 

O' 

z 

* 

ac 

UI 

P-  P- 

z 

—4 

z 

o 

O' 

UI 

• 

uj  a. 

z 

• 

X  >- 

< 

< 

o 

I 

O' 

Q 

• 

ac 

o 

ouu 

z  z 

P- 

—4 

p- 

o 

g 

o 

• 

3  UJ 

Ui 

— « 

Ui 

o 

< 

O' 

-J 

CO 

Ui  > 

p- 

OJ 

1 

M 

*-* 

O' 

UJ 

O' 

< 

p- 

*/»  <t 

a 

Ui  Ui 

a 

• 

O 

• 

• 

ac 

z 

UJ  x 

O 

UJ 

UJ  UJ 

o 

o 

X 

o 

p- 

< 

ac 

CL 

p-  p- 

UJ 

UI 

P- 

g 

p- 

a.  uj 

of 

CM 

< 

<  < 

z 

z 

M 

UI 

Ui 

g 

3 

z  z 

•— « 

X 

a 

z 

u.  <T 

O 

ui 

-J 

*a> 

—I 

CO  CM 

o 

Ui 

o 

O  ac 

u. 

■ 

UJ 

o  a 

co 

r»  o 

3 

o 

o 

o 

p- 

3 

OC  OC 

UJ 

1 

UI 

1 

UI 

CL 

> 

o 

-J 

o  o 

a 

UJ 

Q 

UJ  ui 

P- 

UJ 

UJ 

UJ 

Q 

P-  z 

z 

—t 

<1 

o  o 

o 

O 

O 

O  O 

< 

o 

X 

o 

UJ 

—•  UJ 

> 

g  o 

z 

o 

z 

o  o 

O' 

o 

o 

p- 

o  UJ 

X 

l 

o 

1 

o  o 

o 

O' 

a. 

o 

3 

O  X 

O 

CL 

U_  u 

o 

o 

o 

o  o 

o 

O' 

o 

a. 

— J  p- 

-1 

Ui 

UJ 

O  O 

p- 

o 

p- 

o  o 

CO 

O' 

UI 

o 

X 

aj  uj 

-1 

w 

p- 

1 

o 

1 

o  >r 

VT 

o> 

ac 

o 

o 

>  CD 

o 

z 

UI 

oc  ac 

z 

o 

z 

O  IO 

< 

O' 

3 

g 

u. 

O 

X 

UJ  OJ 

• 

•  • 

• 

CO 

• 

• 

UJ  UJ 

M 

< 

CO  CO 

o 

o 

C3 

o  o 

< 

o 

CO 

o 

• 

O  -J 

UJ 

P> 

oc 

X  X 

•-* 

UJ 

UJ 

* 

<  O 

X 

a. 

< 

3  ZD 

oc 

oc 

ac 

ac 

• 

ac  z 

P- 

O 

a. 

z  z 

o 

O 

< 

a. 

• 

p-  <. 

151 


THE  FOLLOWING  FORCE  CROSS  POWER  SPECTRAL  DENSITY  MATRICES  WERE  GENERATED 


FREOUENCY-  0.3879000E  02(RAO/SEC> 
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DEFLECTION  CO-POWER,  FREQUENCY  -  0.3879000E  02  (RADVSEC! 
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STRESS  CROSS-PSD  MATRICES 
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A  user’s  guide  is  presented  for  a  computer  program  developed  to  aid  in  the  design  of  sonic- 
fatigue-resistant  aircraft  structure.  The  program  employs  matrix  methods  to  calculate 
statistical  measurements  of  response  (deflection  and  stress)  for  complex  structure 
subjected  to  pressure  loads  random  in  both  time  and  space.  The  program  is  in  two  phases. 
Finite-element  methods  are  used  in  the  first  phase  to  determine  structural  characteristics 
such  as  flexibility,  natural  frequencies,  and  modes  of  vibration.  In  the  second  phase,  a 
cross-power  spectral  density  loading  function,  is  generated  and  combined  with  structural 
characteristics  to  compute  response.  Either  cross  power  spectral  density  or  joint  statistical 
moments,  including  second  spectral  moments  useful  in  fatigue  analysis,  can  be  computed 
for  response.  The  loading  function  models  a  decayed  progressive  wave  typical  of  laboratory 
noise  sources.  Different  loading  functions  can  be  supplied  by  the  user,  because  the  program 
is  constructed  in  modular  form.  The  program  was  written  for  the  IBM  7094  computer 
primarily  in  FORTRAN  IV  language  with  a  MAP  language  matrix  manipulation  module. 
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